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ABSTRACT 

We present high-quality, wide spectral coverage long-slit optical spectra for 12 powerful 
radio sources at low and intermediate redshifts (z < 0.7) that show evidence for a 
substantial UV excess. These data were taken using the WHT and VLT telescopes with 
the aim of determining the detailed properties of the young stellar populations (YSPs) 
in the host galaxies as part of a larger project to investigate evolutionary scenarios for 
the AGN host galaxies. The results of our spectral synthesis model fits to the spectra 
highlight the importance of taking into account AGN-related components (emission 
lines, nebular continuum, scattered light) and reddening of the stellar populations in 
studies of this type. It is also clear that careful examination of the fits to the spectra, 
as well consideration of auxilary polarimetric and imaging data, are required to avoid 
degeneracies in the model solutions. In 3 out of the 12 sources in our sample we find 
evidence for broad permitted line components, and a combination of AGN-related 
continuum components and an old (12.5 Gyr) stellar population provides an adequate 
fit to the data. However, for the remaining 9 sources we find strong evidence for YSPs. 
In contrast to some recent studies that suggest relatively old post-starburst ages for 
the YSPs in radio galaxies (0.3 - 2.5 Gyr), we deduce a wide range of ages for the YSPs 
in our sample objects (0.02 - 1.5 Gyr), with ^--^50% of the sample showing evidence 
for young YSP ages (<0.1 Gyr) in their nuclear regions. The nuclear YSPs are often 
significantly reddened (0.2 < E{B — V) < 1.4) and make up a substantial fraction 
(~1 - 35%) of the total stellar mass in the regions sampled by the spectroscopic slits. 
Moreover, in all the cases in which we have sufficient spatial resolution we find that 
the UV excess is extended across the full measureable extent of the galaxy (typically 
5-30 kpc), suggesting galaxy- wide starbursts. The implications for photometric and 
spectroscopic studies of active galaxies are discussed. 

Key words: galaxies: active, galaxies: radio galaxies, galaxies: starburst, galaxies: 
stellar content 



1 INTRODUCTION 

Major merger events are often invoked to explain the origin 
of nuclear activity in galaxies. Indeed, there is much obser- 
vational evidence to suggest that, in a significant fraction 
of the local radio source population, the activity has been 
triggered by merger events involving two or more galax- 
ies where at least one is gas rich. Despite the fact that 
the overwhelming majority of radio galaxies are classified 
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as elliptical or SO types, deep optical imaging studies re- 
veal double nuclei, tidal tails, arcs of emission and dis- 
torted isophotes characteristic of mergers in ~50% of nearby 
radio galaxies with strong em ission lines (Hcckman ot al.l 
1 19861 : ISmith fc HeckmanI 1 19891 ). In addition, kinematical 
studies of the emission line gas (e.g. iTadhunter et al.lll989l : 
iBaum et al] 1 19921 ) further support the idea that the early- 
type hosts are interacting and/or merging with companion 
galaxies. However, there is a limit to what can be deduced 
from the imaging and kinematical studies alone, especially in 
terms of gauging the type of merger (major or minor? dry 
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or gas-rich?), the order of events in the triggering of the 
activity, and the relationship between radio galaxies and 
other types of merging systems such as ultraluminous in- 
frared galaxies (ULIRG). 

Given that gas-rich mergers can also trigger starbursts, 
studies of the stellar populations in the early-type host 
galaxies of radio sources have the potential to provide key 
information about the natures and timescales of the trigger- 
ing events. A major challenge in such studies is to identify 
the young stellar populations in the presence of activity- 
related components such as nebular continuum, direct or 
scattered AGN light and emission lines, all of which are 
likely to b e significant in galaxies hosting powerful, quasar- 
like AGN (jTadhunter et al.ll2002l ). Therefore, although the 
broad-band colours of many radio galaxies show UV excesses 
compared to n ormal, passively evolving elliptical galaxies 
l|Lillv fc Longa ir 1984; Smith & Hcckman 198S), such ex- 
cesses may be related more to the AGN activity than any 
merger induced star formation. Nonetheless, by taking par- 
ticular care to account for all the activity related compo- 
nents, it has proved possible to use spectral synthesis mod- 
elling techniques to reveal the presence of young stellar 
popu lations (YSPs ) in radio galaxie s (e.g. IXadhunter et al.l 
Il996h . In this way IXadhunter et al.1 ()2002l ) found evidence 
for YSPs in 30 - 50% of a complete sam ple of southern 2Jy 
radio galaxies at intermediate redshifts. IWills et al.l (|2002l ) 
deduced the presence of YSPs in 30% of their complete 
sample of nea rby 3C radio galax ies with FRII radio mor- 
phologies, and IWills eral] (120041 ) detected YSPs in 25% of 
a small sample of lower radio power FRI radio source s from 
the 2Jy sample. More recentlv. iRaimann et al ](2005) found 
a higher rate of incidence of YSPs in a sample of 20 nearby 
radio galaxies (~100%). However, the latter authors did not 
explicitly account for all the AGN-related continuum and 
emission line components, and assumed a priori that the 
optical spectra of their radio galaxies could be modelled in 
terms of 5 stellar components of varying age. 

Although all the previous studies revealed the presence 
of young stellar populations in a significant fraction of radio 
galaxy hosts, the data were not always of sufficiently high 
quality in terms of spectral range, spectral resolution and 
S/N to allow the detailed properties of the YSPs to be de- 
duced, and thereby learn about the nature of the triggering 
merger events. We have therefore begun a major programme 
to obtain the requisite data and investigate the stellar popu- 
lations in all nearby radio galaxies already suspected of hav- 
ing YSPs on the basis of previous work, aiming to analyse 
the data for all the objects in a uniform way. Some early re- 
sults from this progra mme for 5 nea r by radio galaxi e s were 
pres ented in .Tadhuiiter et al.l (|2005h : lEmonts etall (120061 ) 
and I Holt et aL ( 20061 ) . The latter papers also give a full ac- 
count of the background to the project. In this paper we 
present the results for a further 12 radio galaxies at low and 
intermediate redshifts [z < 0.7j3, extending detailed studies 
of the YSPs in radio galaxies up to higher redshifts than the 
previously published work. The implications of the results 
of the programme as a whole for our understanding of the 
evolution of the host galaxies, triggering of the activity and 
the heating of the far-IR emitting dust in AGN are discussed 



Ho = 75 km s ^ Mpc ^, go = 0.0 assumed throughout. 



in separate papers (jTadhunter et al.ll2007l , Tadhunter et al. 
in prep). 



2 THE SAMPLE 

The sample considered in this paper comprises 12 powerful 
low and intermediate redshift [z < 0.7) radio galaxies. All 
have been selected on the basis of having a strong UV excess 
that is difficult to explain solely in terms of AGN-related 
components (see the discussion in Tadhunter et al. 2002) . In 
essence, all sources in this sample are suspected of having a 
large contribution of continuum light from a young stellar 
population (YSP). Therefore they are suitable for detailed 
studies of the young stellar populations. 

The sample spans both the northern and southern skies. 
The northern objects originate from the 3C catalogue whilst 
the southern objects derive from the well-studied 2Jy sample 
of radio galaxies (^c.g. Tadhunter ct al.''l99A lMorganti et al.l 
[1993; Tadhun ter et al, 2002 ). Table [1] lists the properties of 
the sample. 

Although many of the objects in our sampl e have been 
previously modelled ( Tadhunter et al.i ri996: Ar etxaga et al 

^13, 



20041 : IRaimann et al 



200ll : lRobinsonl200ll : [wills et al.ll20o' 
20051 ). for most we have new data with improved spectral res- 
olution, spectral coverage and S/N. We also fit these data 
with the newer, higher resolution spectral synthesis models 
of Bruzual & Chariot (2003), taking full account of redden- 
ing of the YSPs and potential AGN contributions to the 
optical/UV continuum. 



3 OBSERVATIONS AND DATA REDUCTION 

The northern sample was observed using the ISIS dual-beam 
spectrograph on the William Herschel Telescope (WHT) on 
La Palma in several runs between 1996 and 2005. The south- 
ern sample was observed during a single run in 2003 using 
the FORS2 spectrograph on the ESQ Very Large Telescope 
(VLT) on Cerro Paranal, Chile. Table [5] summarises the ob- 
servations and the various instrumental setups. Note that, 
in the case of the highest redshift source in our sample (also 
the faintest) — PKS 0409-75 — we nodded the source along 
the slit in order to improve the sky subtraction at the red 
end of the spectrum. Al though for 3C 3 2 1 we use the same 
data set as that used in iRobinson et al.l (|2000l ). for this pa- 
per, we have improved the flux calibration and re-modelled 
the SEDs using the latest models, for the first time taking 
into account reddening of the YSP. 

This project focusses on accurate modelling of the (of- 
ten faint) continuum emission. Hence, all observations were 
made with a 1.3-1.5 arcsec slit in order to include as much of 
the light from the galaxy as possible, whilst retaining suffi- 
cient wavelength resolution to allow detection and modelling 
of stellar absorption features. To reduce the effects of differ- 
ential refraction, all exposures were taken at low airmass (sec 
z < 1.1), with the slit aligned along the parallactic angle, 
and/or (in the case of the VLT observations) using a linear 
atmospheric dispersion compensator (LADC) that was reset 
between exposures. The slit PAs are listed in Table [2] (for 
reference, the radio axis PAs are presented in Table [T|. 
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Table 1. Properties of the sample, (a) radio source; (b) redshift; (c) Galactic reddening from Schlegel et al. 
(1998) (E(B-V)); (d) radio luminosity at 5GHz calculated from the flux densities in Wall & Peacock (1985) 
except 3C 285 (from Kellerman et al. 1969); (e) radio spectral index for oI'^qIj^ where oc u~°', (/) radio 

morphology and (g) optical spectral classification (see text fo r references). 

redshifts: * Holt (2 005); f van Breugel fc D ey (1993); o Ide Grijp et al.l lll992l) ; o iHill et al.l (1996); t> 
ISmith et al * iTadhunter et al ] l|l993l ): < Ide Vaucouleurs et al.l l|l99lh . References to the other data 

can be found in the text. 



z Galactic Radio a Radio Spectral 

reddening luminosity (5 GHz) morphology class 

(E(B-V)) log P (W Hz-l) 

(«) (b) (c) id) (e) if) (g) 



3C 218 


0.05488> 


0.042 


25.91 


0.9 


FRI/FRII, symmetric 


NLRG 


3C 236 


O.lOlo 


0.011 


25.44 


0.7 


FRH, double-double 


NLRG 


3C 285 


0.079t 


0.017 


25.00 


1.27 


FRH 


NLRG 


3C 321 


0.096O 


0.044 


25.33 


1.14 


FRH 


NLRG 


3C 381 


0.161t 


0.053 


25.87 


0.96 


FRH 


BLRG 


3C 433 


O.lOOf 


0.044 


25.94 


1.07 


FRI/FRH, x-shaped 


NLRG 


PKS 0023-26 


0.322,t 


0.014 


27.33 


0.7 


CSS 


NLRG 


PKS 0039-44 


0.346* 


0.007 


26.57 


0.93 


FRH 


NLRG 


PKS 0409-75 


0.693* 


0.078 


27.82 


0.86 


FRH 


NLRG 


PKS 1932-46 


0.231* 


0.054 


26.65 


1.03 


FRH 


BLRG 


PKS 2135-209 


0.636* 


0.033 


27.47 


0.82 


CSS 


BLRG 


NGC 612 


0.030<] 


0.020 


24.85 


0.51 


FRI/FRH, 'hybrid' 


NLRG 



3.1 Data reduction 



4 RESULTS 



The data were reduced in the usual way (bias subtraction, 
flat fielding, cosmic ray removal, wavelength calibration, flux 
calibration) using the standard packages in iraf. The two- 
dimensional spectra were also corrected for spatial distor- 
tions of the CCD. The final wavelength calibration accura- 
cies, calculated using the standard error on the mean devi- 
ation of the night sky emission line wavelengths from pub- 
lished values (Osterbrock et al.lll99d ). ranged between 0.05- 
O.5A. The spectral resolution, calculated using the widths 
of the night sky emission lines, ranged between 4-loA. Both 
are dependent on the instrumental setup and the wavelength 
range observed. 

Comparison of several spectrophotometric standard 
stars observed with a wide slit (5 arcsec) throughout each 
run gave relative flux calibrations accurate to ±5 per cent. 
This accuracy was confirmed by good matching in the flux 
between the red and blue spectra. For the VLT observa- 
tions, and for 3C 236, further observations of standard stars 
close in position and time to each target observation with a 
narrow slit, matched to the slit width used to observe the ob- 
jects, were used to correct for atmosperic absorption features 
(e.g. A and B bands at ~7600 and ~6800 A respectively). 
The spectra of all sources were corre cted for Galac t ic ex - 
tinction using the E(B -V) values from lSchlegel et al.l (|l998l ) 
and the lSeatonI (|l979l ) extinction law. 

The spectra were extracted and analysed using the 
STARLINK packages figaro and dipso, and confix, a cus- 
tomised IDL minimum fltting programme (|Robinsonl200ll : 
ITadhunter et al.ll2005h . 



The continuum spectra for all sources were analysed using 
three main steps: first, the degree of the UV excess and 
its spatial distribution were assessed using the 4000A break; 
second, the overall spectral energy distribution (SED) shape 
was modelled using a combination of nebular, stellar and 
power-law continua (which could represent an AGN-like or 
very young stellar population continuum); and finally, mak- 
ing detailed comparisons between the better fitting SED 
models and the data using stellar absorption lines (e.g. 
Balmer lines, Call K). The individual steps are described 
in detail below. 



4.1 The 4OOOA break 

A simple, commonly used diagnostic tool for investigat- 
ing the magnitude of the UV excess in a galaxy spectrum 
is to compare th e continuum flux on either side of the 
4000 A break (e.g. ITadhunter et al.|[20o3 . |2005| ). Hence, for 
each galaxy, continuum slices on either side of 4000A were 
extracted to estimate the D(4000) ratio - a measure of 
the 4000 A brea k amplitude - following the definition in 
ITadhunter erall l|2002l . l2005h : 



l,,r„ F\d\ 
J3750 -P^^d^ 



(1) 



This form of the 1)4000 ratio is be tter suited to active galax - 
ies than the original form (see e.g. iBruzual" fc Charlot|[l993l ) 
as it is not contaminated by bright emission lines (e.g. [Ne 
III]A3869). The L>(4000) ratio is also insensitive to both 
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Figure 1. D4000 ratio for all galaxies in the sample. The solid and dashed lines (profiles) represent the integrated continuum flux in the 
ranges 3750-3850A and 4150-4250A respectively. The D(4000) ratio is overplotted (faint solid line) along with the expected ratios for a 
pure 12.5 Gyr elliptical galaxy (bold solid horizontal line) and diluted with a 50 per cent featureless continuum in the 3800A bin (dashed 
horizontal line) and with a 75 per cent featureless continuum (dotted horizontal line). The ratio expected for a subsolar metallicity (Z = 
0.4 Zq) 12.5 Gyr population is also shown (dot-dashed horizontal line). The locations of all apertures used in the analysis are also marked. 
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Table 2. Summary of the observations. Note, the wavelength ranges given are the useful wavelength ranges of the spectra. Further, 
where the end/start wavelengths of the blue/red frames are equal, these are where the 'trim' was made for the fitting - in the 
majority of sources there was considerable overlap between the blue and red frames, of order a few hundred angstroms. 



Date 


Object 


Setup 


Exposurc 




Slit 


Rest A 


Seeing 






(arm / CCD / grating/filter) 




PA 


width 


range 












(°) 


(arcsoc) 


(A) 


(arcsec) 


WHT/ISIS 
















13/01/2004 


3C 218 


B /EEV12 /R300B /- 


3*1200 


10 


1.27 


3000-5950 


~1 






R/MARCONI2 /R316R/GG495 


3*1200 


10 


1.27 


5950-8000 




12/05/2001 


3C 236 


B /EEV12 /R300B /- 


2*1200 


50 


1.3 


3040-5640 


0.7 






R /TRK4/R316R /CG495 


2*1200 


50 


1.3 


5640-6880 


0.7 


06/06/2005 


3C 285 


B /EEV12 /R300B /- 


3*1200 


209 


1.23 


3100-4848 


0.6 






R/MARCONI2/R316R/- 


3*1200 


209 


1.23 


4848-6994 


0.6 


22/07/1996 


3C 321 


B /LORALl /R158B /- 


2*1500 


130 


1.37 


3309-5700 


<1.0 






R/TEK5 /R158R/GG495 


2*1500 


130 


1.37 


5883-8100 


<1.0 


22/07/1996 


3C 381 


B /LORALl /R158B /- 


2*900 


20 


1.34 


3123-5560 


0.8-1.0 






R/TEK5 /R158R/GG495 


1*900,1*1300 


20 


1.34 


5560-8464 


0.8-1.0 


31/07/2005 


3C 433 


B/EEV12/R300B/- 


3*1200 


32 


1.54 


2800-4900 


1.1-1.7 






R/MARCONI2/R316R/GG495 


3*1200 


32 


1.54 


4900-7433 


1.1-1.7 


VLT/FORS2 














24/09/2003 


PKS 0023-26 


B/MIT/G600B/- 


3*1000 


145 


1.3 


2850-3800 


0.7-0.8 






R/MIT/G600RI/GG435-I-81 


3*1000 


145 


1.3 


3800-6100 


0.7-0.8 


25/09/2003 


PKS 0039-44 


B/MIT/G600B/- 


3*1000 


72 


1.3 


2900-4200 


0.7-1.0 






R/MIT/G600RI/GG435+81 


3*1000 


72 


1.3 


4200-6132 


0.7-1.0 


24/09/2003 


PKS 0409-75 


R/MIT/G600RI/GG435-f81 


12*600 


-30 


1.3 


2921-4875 


1.0-1.2 






R/MIT/G600Z/OG590+32 


6*600 


-30 


1.3 


4875-5800 


1.0-1.2 


24/09/2003 


PKS 1932-46 


B/MIT/G600B/- 


3*900 


9 


1.3 


3000-4900 


0.7-1.0 






R/MIT/G600RI/GG435-I-81 


3*900 


9 


1.3 


4900-6700 


0.7-1.0 


25/09/2003 


PKS 2135-20 


B/MIT/G600B/- 


3*900 


52 


1.3 


2250-3600 


1.3-1.6 






R/MIT/G600RI/GG435-I-81 


3*1200 


52 


1.3 


3600-5000 


1.3-1.6 


24/09/2003 


NGC 612 


B/MIT/G600B/- 


3*900 


-13 


1.3 


3500-5855 


0.6-0.9 



reddening and redshift (i-f-corr ection) effects, unlike bro ad- 
band colour comparisons, (e.g. ISmith fc HeckmanlllQigl '). 

Figure [1] shows the spatial distribution of the con- 
tinuum flux in the two regions around 4000A for all 
sources in the sample along with the derived _D(4000) ra- 
tio. For comparison, expected -D(4000) ratios in the 3750- 
3850A bin are overplotted . These ratios were derived us- 
ing the iBruzual fc CharlotI (|2003l ) spectral synthesis mod- 
els (hereafter BC03) for a 12.5 Gyr old stellar population, 
formed in an instantaneous burst with solar metallicity and 
a Salpeter initial mass function (IMF). The various com- 
parison 73(4000) ratios comprise: a) a 12.5 Gyr elliptical 
with solar metallicity, b) a 12.5 Gyr elliptical with sub-solar 
metallicity (Z — 0.4 Zq), c) a 12.5 Gyr elliptical diluted with 
a fiat spectrum such that the 12.5 Gyr population accounts 
for i) 50% and ii) 25% of the continuum flux. 

All sources in the sample show a significant UV ex- 
cess compared to a passively evolving old (12.5 Gyr), solar 
metallicity elliptical galaxy. For all galaxies in the sample, 
the D(4000) ratio suggests that only 25-50% of the flux in 
the UV originates from the 12.5 Gyr stellar population. In- 
terestingly, for those sources that are spatially resolved (3C 
236, 3C 218, 3C 321, 3C 285, 3C 433, PKS 0023-26, NGC 
612), the UV excess is also spatially extended, covering the 
entire spatial extent of the galaxy, suggesting that, whatever 
causes the UV excess (AGN continuum and/or a YSP) is 
important galaxy wide, and is not confined to the nuclear 



regions. Further, some objects show interesting trends in 
their UV excess. In 3C 218, the UV excess increases sym- 
metrically towards the centre of the galaxy. In constrast, the 
UV excess in NGC 612 decreases symmetrically towards the 
nucleus. It is also notable that, in the case of the double 
nucleus systems 3C 321 and 3C 433, the UV excesses ex- 
tend across both nuclei as well the extended haloes of the 
systems. 

The spatial continuum profiles in Figure [T] have been 
used to identify regions of interest in the 2-dimensional spec- 
tra for extraction and continuum modelling (see below) . For 
the majority of objects, only a central nuclear aperture was 
extracted. Notable exceptions are NGC 612, 3C 433 and 3C 
321. All apertures are marked on Figure [T] Our complete 
spectra cover a wide wavelength range, vital for accurate 
modelling of the SEDs, and the observed rest wavelength 
ranges for each object are summarised in Tabled 



4.2 The nebular continuum 

Prior to continuum modelling, a nebular continuum was gen- 
erated for each aperture to be analysed over the entire wave- 
length range. This comprises the blended higher Balmer se- 
ries (> H8) and a theoretical nebular continuum (gener- 
ated using the starlink package dipso, which includes the 
effects of free-free emission, free-bound recombination and 
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two-photon continua) following iDickson et al.l (|l995l ). Both 
components of the nebular continuum were normalised to 
the H/3 flux in each emission line component. Note, H/3 was 
modelled using the minimum number of Gaussian compo- 
nents required to give a physically viable good fit and so 
for many apertures, in particular the nuclear aperture, the 
generated nebul ar continuum inclu ded multiple components, 
as discussed in iHolt et aL I (12003') for the complex nuclear 
aperture in PKS 1345-1-12. The spectra of several of the ob- 
jects are subject to significant reddening. Hence, the nebu- 
lar continuum was generated using the de-reddened H/3 line 
fluxes and then re-reddened before being subtracted from 
the galaxy spectrum. The degree of internal reddening was 
estimated using the Balmer line decrements (e.g. Ha/H/J, 
H7/H/3) and estimated separately for the different kinematic 
components of the lines in each aperture. 

It should be noted that, when generating the nebu- 
lar continuum, one must carefully consider the potential 
pitfalls. In sources where a significant fraction of the op- 
tical/UV light originates from a young stellar component, 
the Balmer lines (particularly H/3, H7) may be affected by 
underlying stellar absorption lines, and the flux therefore 
under-estimated. To check on the likely impact of the stel- 
lar ab sorption lines, we have followed the technique of iHoltl 
l|2005l ) and subtracted off the best fltting continuum model 
and then re-measured the emission lines. 

Considering these issues, the checks show that the spec- 
tra fall into three categories: i) apertures in which the origi- 
nal and re-calculated nebular continua based on Ha/H/3 are 
similar (within a few percent in the 3450-3550A bin; e.g. 
PKS 0039-44) - we flnd that, due to this effect, H/3 is typi- 
cally under-estimated by < 10 per cent. Further, in apertures 
where reddening is important, a 10 per cent correction on the 
H/3 flux, gives a variation in the E(B-V) value of only <0.1 
for Hq/H/3 < 10. For a discussion regarding the extinction 
law used, see Section [5.21 ii) apertures in which the nebular 
continuum fraction is so small (< few per cent) that even 
large differences in nebular reddening would not significantly 
affect the results (e.g. 3C 218, 3C 236, NGC 612) and in) 
apertures in which the reddening estimate was based solely 
on weaker lines (e.g. H7, H/3) due to Ha lying outside of the 
observed spectral range (e.g. PKS 0023-26, PKS 2135-20, 
PKS 1932-46). For case iii), we are less confident of the red- 
dening estimate due to uncertainties in measuring H7 as it is 
a weak, blended line which is relatively more affected by un- 
derlying stellar absorption features than Ha or H/3. In these 
apertures, we h ave followe d the example of PKS 1345-1-12 
as discussed in iHoh et al.l |2OO30 which allows us to de- 
termine the maximum nebular continuum fraction, i.e. that 
which does not produce an unphysical step-function in the 
data at the location of the Balmer break when subtracted. 
For these apertures, we have modelled the continuum twice, 
with both the maximum and minimum (i.e. 0%) nebular 
continuum fraction (this technique was also used for PKS 
1549-79; see H06). Reassuringly, we find that the modelling 
results are similar for both cases of nebular continuum sub- 



^ As the technique for correct ing for the n ebula r continuum used 
here is discussed in detail in iHoIt et al.l , we do not dis- 

cuss the technique or show examples of the generated nebular 
continua. 



traction (minimum and maximum). Assuming a zero nebu- 
lar continuum fraction gives a slightly larger UV continuum 
fiux and so would potentially shift the Xrod minimum to 
slightly younger ages. For apertures in which the continuum 
is dominated by a power-law component (e.g. PKS 1932-46) 
and/or the YSP age is poorly constrained (e.g. PKS 2135- 
20), the modelling results based on both the minimum and 
maximum nebular continua are indistinguishable i.e. the ef- 
fect on our results due to errors on the subtracted nebular 
continuum are insignificant. In the nuclear aperture of PKS 
0023-26, the YSP contribution is large (50-60%) and has a 
tightly constrained minimum (see Figure (J)). For this case, 
subtracting off the maximum nebular continuum (22.5%) 
gives a YSP age range of 0.03-0.05 Gyr whilst a zero neb- 
ular continuum model gives a YSP age range of 0.01-0.03 
Gyr. We are therefore confident that the uncertainties in 
the nebular continuum subtraction do not significantly af- 
fect our results. 

The contribution of the nebular continuum (reddened 
where necessary) to the fiux in the 3540-3640A bin in each 
aperture ranges from <1% (NGC 612, 3C 381, mid aperture 
in 3C 321) to ~25% (extended aperture of 3C 321, PKS 
1932-46) and is summarised in Table [S] 

For the objects with published nebular continuum mea- 
surements, the nebular continuum fraction previously es- 
timated is typically higher than the estimates here. For 
the sources in the 2Jy sample (PKS 0023-26, PKS 0039- 
44, PKS 0409-75, PK S 1932-46 and PKS 2135-20 - see 
iTadhunter et al.]|2002l). the differenc es are likely to be due to 
the fact that Tadhunter et al.| (|2002l ') made no correction for 
reddening, stating that these nebular continuum fractions 
were upper limits. In this paper, we have followe d the care- 
ful re ddening correction technique described in iHolt et al] 
(l2003l ) and generally derive lower nebular continuum frac- 
tions. 



4.3 Continuum modelling 

After subtracting the nebular continuum, the continuum 
SEDs were modelled using CONFIT, a purpose writte n idl 
code using a minimum technique l|Robinsonll200ll ). The 
program allows up to three continuum components to be 
combined - up to two different stellar populations, with 
or without a power-law component. Note that the power- 
law could represent an AGN/scattered light component or 
a very young stellar population. The stellar templates used 
in the modelling were taken from the 2003 version of the 
spectral synthesis results of Bruzual & Chariot (2003; here- 
after BC03). The stellar population ages used include old, 
evolved stellar populations (OSP: 12.5 Gyr and, in some 
cases, 7.0 Gyr) and various young stellar populations (YSP: 
0.01-5.0 Gyr in varying steps). All stellar population models 
assume a Salpeter initial mass function (IMF), solar metal- 
licity, instantaneous starburst. The YSP spectra were red- 
dened ov er the range 0.0 < E(B-V)< 1.6 in steps of 0.1 
using the lSeatonI l| 19791 ) extinction law (see Section [5. 2|) . Fi- 
nally, when a power-law component was included (of the 
form Fa oc A"), the spectral index of the power-law was al- 
lowed to vary between -15.0 and 15.0 (to ensure reasonable 
calculation times). 

Various models were tried including: 
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UV excess? 

(from D4000 ) 



Yes 



No 



Unlikely to have 
large contribution 
from YSP 




SED modelling 
absolutely requires 
YSP to adequately 
fit SED and 
absorption lines 



\ 




Likely to have a 
significant YSP 



What is the age/ 
reddening/mass of 
the YSP component? 





SED adequately 
fit with power-law 
andOSP 








V 



If fitted power-law is red 
it is unUkely to represent 
a scattered AGN - in this 
case, failure to detect HST 
point source (expected for 
direct AGN Ught), would 
strongly favour a reddened 
YSP 



Large polarisation 
inUV? 



No 



Broad permitted 
Unes detected? 
(e.g. Ha, HP, Mg II) 



Yes 



Likely to be dominated 
by scattered AGN or 
non-thermal continuum 
source 



No 



Yes 



Point source detected 
in optical HST image? 



UV excess Ukely to 
be dominated by 
direct AGN Ught 
(weakBLRG/QSO) 



No 



Yes 



Strong emission lines 
with large EW? 



No 



Yes 



UV excess may be 
dominated by direct AGN 
light (but low level 
because no broad lines 
or a large UV polarisation 
is detected) 



UV excess likely to 
be dominated by 
YSP with very yoimg 
age (+ reddening?) 



UV excess may be 
dominated by a very 
yoimge age YSP, but 
residual low-level AGN 
(direct or scattered) 
also possible 



Figure 2. The 'decision tree' used to interpret the modelling results. Abbreviations arc as follows: SED: spectral energy distribution; YSP: 
young stellar population; OSP: old stellar population; EQW: equivalent width; BLRG: broad line radio galaxy; QSO: quasar; 1)4000: ratio 
of continuum fluxes at the 4000A break. 



• an OSP (12.5 Gyr) with zero reddening. 

• an OSP (12.5 Gyr, zero reddening) with a power-law 
component. 

• an OSP (12.5 Gyr, zero reddening) with a YSP (0.01 < 
tysp < 5.0 Gyr, 0.0 < E(B-V)< 1.6). 

• an OSP (12.5 Gyr, zero reddening) with a YSP (0.01 



< tysp < 5.0 Gyr, 0.0 < E(B-V)< 1.6) and a power-law 
component. 

• for the higher rcdshift sources (e.g. PKS 2135-20, PKS 
0039-44, PKS 0409-75 and PKS 0023-26), 'younger' old stel- 
lar populations were also tried for all four combinations 
above including 7.0 Gyr, 9.0 Gyr and 10.0 Gyr because, at 
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the redshifts of these objects, current cosmological models 
show the universe is younger than 12.5 Gyr. 

Our general philosophy is to model the data with the 
minimum number of stellar and/or power-law components 
that provide an adequate fit to the spectra. Given that 
the galaxies we are modelling are all early-type in terms 
of their optical morphologies, we start by assuming that all 
galaxies contain an old (7 - 12.5 Gyr) stellar population 
(OSP) at some level, and that any young stellar popula- 
tions were formed in a single instantaneous burst as a result 
of a merger, and hence have a particular age. In our models 
very young, reddened YSPs (<10 Myr) may be simulated 
by adding a power-law component, provided that they can 
be distinguished from AGN-related components (see below) . 
Clearly more components would be required to model more 
complex star formation histories, but the greater the num- 
bers of components, the greater the likelihood of degenera- 
cies in the solutions. 

One of the major ambiguities in our model solutions 
is that we cannot readily distinguish between the case of 
a highly reddened very young YSP (<10 Myr) plus an 
OSP, and an AGN-related power-law plus an OSP, based 
on the SED modelling alone. This is because the spectra of 
the very young YSPs are relatively featureless in terms of 
the strengths of the a bsorption lines and the ampli t ude o f 
the Balmer break (e.g. lGonzalez Delgado et al.ll 19991 . l2005h . 
Therefore it is necessary to use auxilary information, such as 
UV polarisation, detection of point sources in HST images, 
detection of broad permitted lines, and emission line lumi- 
nosity to distinguish between these cases. For example, if 
the UV excess is predominantly caused by a scattered AGN 
component, revealed as a power-law in the SED fits, we ex- 
pect to measure significant UV polarisation and a relatively 
fiat or blue power- law slope (a < 0.5j3- On the other hand 
we expect a direct AGN component (from a low-luminosity 
or partially obscured quasar) to be associated with the de- 
tection of broad permitted lines and a point source in HST 
images. Figure [2] presents a flow diagram of the decision 
making process involved in interpreting the modelling re- 
sults for such cases. 

For the modelling, continuum bins were selected to 
avoid emission lines and, in the case of 3C 285, 3C 433, 3C 
218/Hydra A, 3G 321 and 3G 381, the uncorrected atmo- 
spheric absorption bands. Typically, 40-50 bins were chosen 
for each aperture to ensure good sampling over a large spec- 
tral range as evenly distributed in wavelength as possible. 
In addition, a normalising bin in the blue part of the spec- 
trum, common to all objects, was chosen with rest wave- 
length 4720-4820 A. The confit program then generated 
the models by scaling the different components to the flux 
in the normalising bin so the total model flux was <125 per 
cent of the observed flujfl, flnding the model combination 
with the minimum reduced (Xrcd)- For the x?cd determi- 



Such blue power-law slopes have been meas ured for all cases 
with significant scattered AGN components llTad hunter et al.l 
1 19961 : iTran et al.lll998l : ICohen et ai]|l999l : IVernet et al.ll200ll) and 

strong emission lines — the objects with the strongest scattered 
AGN components tend to bo those with the most luminous emis- 
sion lines (see fradhuntcr ct al. 2002) 

* The maximum allowed model flux is greater than 100 per cent 



nation, we assume a relative flux calibration error of ±5 per 
cent. 

Gombining the best fitting model for each combination 
of OSP with reddened or unreddened YSPs of various ages 
(both with and without a power-law component) forms a 
surface in xLd for YSP age versus E(B-V) from which con- 
tour plots can be generated. Examples of the contour plots 
of the x?od space for each aperture requiring a YSP compo- 
nent are shown in Figure [3] Using these plots, we identified 
minimcffl in the x?od space to be further investigated us- 
ing detailed comparisons of the stellar absorption features 
known to be sensitive to age. The Gall KA3968, the G-band 
(4306A) and the Mg lb band (5173A) absorption features 
are particularly important for such comparison because they 
are not significantly infilled by emission lines. In some cases 
with weak emission lines it was also possible to use the higher 
order Balmer line features for detailed comparisons. Along 
with examples of the overall SED fits, detailed comparisons 
of the region containing the Balmer series and Gall H-fK 
are shown in Figure [l] 

Note that, when discounting models on the basis of de- 
tailed fits to Gall K feature, consideration must be made of 
the potential contribution of intrinsic ISM absorption to the 
to this feature, especially in highly reddened cases. A con- 
tribution from the ISM would strengthen the Gall K line 
in the observed spectrum, but not in the model spectra. 
Therefore models that underpredict the observed strength 
of Gall K may be viable, whilst models which significantly 
overpredict the strength of Gall K can be discounted. This 
technique proved particularly useful in distinguishing be- 
tween two clear minima in the x?od space (e.g. 3G 285). In 
some cases it was necessary to apply smoothing to the mod- 
els to obtain good matching in the absorption lines because 
the models and the data have different effective resolutions. 

Hence, in order to be considered a viable fit to the data, 
the various models had to provide both a good fit to the 
overall SED and the stellar absorption lines. When relying 
on one technique alone, the solutions are of ten degenerate in 
age/reddening space for the YSP (see also iTadhunter et al.l 
[200a . 

4.3.1 Hydra A (3C 218) 

3G 218 is one of the most powerful radio galaxies in the local 
universe, however, despite having a radio luminosity more 
than an order of magnitude larger than the characteristic 
break luminosity between the FR I and FR II classes, its 
morphology is edge darkened — more characteristic of FR I 

sources (e.g. iBaum et al]|l988l ). 

3G 218 is optically identified with Hydra A l|Simkinl 

Il979l ). the dominant cD 2 galaxy in the c entre of a poor 
but X-ray luminous (e.g. iDavid et al]|l990t ) cluster. Hydra 
A has a double nucleus with the second nucleus lying ^ 
7 arcsec SE of the active n ucleus, within the common en- 
velope ([Hansen et al.l [l996l ). Like a significant number of 

of the measured flux to allow for uncertainties in the measured 
flux and models in the bin. 

® Note, whilst the main discriminating factor is the minimum in 
the Xrcd space, all fits with Xred ~t are deemed to be a 'good fit'. 
In addition, fits with Xj;,;j^l-2 are also considered for the detail ed 
comparisons. See detailed discussion in ITadhunter et al.l fcOOSl '). 



3C 218, Hydra A (nucleus): 12.5 Gyr + YSP 



3C 236 (nucleus): 12.5 Gyr + YSP 
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3C 285 (nucleus): 12.5 Gyr + YSP 
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3C 321 (SE nuc): 12.5 Gyr + YSP + power law 
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3C 321 (NW nuc): 12.5 Gyr + YSP + power law 



Age (Gyr) 

3C 321 (mid scaled): 12.5 Gyr + YSP + alpha 





Age (Gyr) 
3C 381 (nuc): 12.5 Gyr + YSP 




Age (Gyr) 

3C 433 (SW nuc): 12.5 Gyr + YSP 



Age (Gyr) 

3C 433 (NE nuc): 12.5 Gyr + YSP 



Age (Gyrl 
3C 433 (mid): 12.5 Gyr + YSP 

















NX" \ 






u\\\\\ ■ \\ 



























--- 1.0. 


































Figure 3. Contour plots of the reduced chi squared space (x^ed) ^" combinations of age (x axis) and reddening (y axis) of the YSP 
for all apertures requiring a YSP componet. All plots assume the [Scaton 19791 reddening law except for the second (labelled) plot for the 
nuclear aperture of NGC 612 which uses the lCalzetti et al]|2000l reddening law. For descriptions of the plots see Section 14.31 In general, 
we show one contour plot for each aperture. For apertures where a power-law component is expected, the plot is for models with an OSP, 
a YSP and a power-law component. For apertures in which a power-law component is not expected (e.g. due to location, the weakness of 
emission lines), the contour plot is for models with an OSP and a YSP only. Where the contours in the region of the minimum are crowded, 
a zoom of this region is shown inside the pl ot. Note, a co ntour plot is not shown for PKS 0039-44. Two plots for the n uclear aperture of 
NGC 612 are presented, the first using the l ISeatonlllOTgl) extinction law and the second using the llCalzetti et al.ll200Cll law. These plots 
illustrate how the different reddening laws may influence the result. The plots clearly show that, whilst the details of the models may vary, 
the differences between the reddening laws do not change the overall result within a certain tolerance. We discuss the different reddening 
laws further in Section FS. 2 1 and Figure [Tol 
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PKS 0023-26 (iiuc): 12.5 Gyr + YSP 



PKS 0409-75 (nucleus): 7.0 Gyr + YSP 



PKS 1932-46 (nucleus): 12.5 + YSP + power law 






NGC 612 (nucleus): 12.5 Gyr + YSP (Seaton 1979) 



Age (Gyr) 

NGC 612 (nuc): 12.5 Gyr + YSP (Calzetti 2000) 



Age (Gyr) 
NGC 612 (A): 12.5 Gyr + YSP 





Age (Gyr) 
NGC 612(B): 12.5 Gyr + YSP 
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NGC 612(C): 12.5 Gyr + YSP 



Age (Gyr) 
NGC 612(D): 12.5 Gyr + YSP 






Age (Gyr) 
NGC 612 (E): 12.5 Gyr+ YSP 



Age (Gyr) 
NGC 612 (F): 12.5 Gyr + YSP 




PSK 2135-20 (nucleus): 12.5 Gyr + YSP + power law 
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Figure 3. Contour plots continued. 



radio galaxy hosts, the optical colours of 3C 218 are blue 
in comparison to normal elliptical galaxies (see also Fig- 
ure Several authors have attributed these colours to 
a massive (10^-10^ M(;)) burst of star formation occurin g 
- 50Myr ago (e.g. iMcNamaral Il995l : [Hansen et aLlll996l ). 
Optical spectroscopy supports this view with the detec- 



tion of a prominent Balmer series in the n uclear sp ectrum 
(.Hansen et al. 1996; Mclnick ot ahj 1 19971 : Arctxag a et al.l 
I2OOII : IWills et al.l 12004) Using a variety of line and break 
indices. lAretxaga et al.l (|2001i ') attribute the blue bulge- 
subtracted continuum to a 0.007-0.04 Gyr old starburst 
dominated by B-type giant/supergiant stars. Spectral syn- 
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thesis mode lling by Wills et al. ( 2004 ) finds results con- 
sistent with lAretxaga et alj ~i 2001 ) although their best-fit 
model suggests the starburst occurred at the earlier end of 
this range, 0.05 Gyr ago. 

In this paper we h a ve re peated the optical spectral 
modelling of IWills et al.l (|2004l ') using new, deeper spectra 
with wider spectral coverage, and usin g the la test gen era- 
tion o f high resolution spectral models of Bruzual fc CharlotI 
l|2003ll . 

As discussed above, and shown in Figure [T] 3C 218 
has a clear UV excess and the prominent Balmer series 
m akes it a strong sta rburst candidate. Indeed, as reported 
by I Wills et"aLl (|2004 ). we find that the SED cannot be ad- 
equately modelled without a YSP component. Hence, we 
have modelled the SED with an OSP and a YSP both with 
and without a power- law component. In Figures |3] and |4] we 
only show the case of an OSP plus a YSP - in common with 
many other cases, adding a power-law gives a similar shape 
to the Xred space with a similar position for the x?od rnin- 
imum although the number of solutions giving a viable fit 
to the overall shape of the SED is dramatically increased. 
Whilst including a power-law component gives many viable 
fits, there is no strong observational evidence for a power- 
law component in 3C 218. Despite harbouring a powerful 
radio source, the optical emission lines are relatively weak 
(we have not corrected for the nebular continuum) and there 
is no strong point source in any of the published images (e.g. 
lMcNamaralll995l - ground based). 

For fits without a power-law component, one clear min- 
imum is observed at 0.05 Gyr with significant reddening 
(E(B-V) = 0.4; Xrod = 1-2) accounting for 50% of the con- 
tinuum flux. In fact, only YSPs with ages 0.05-0.1 Gyr with 
similar reddening provide acceptable fits to the SED mod- 
elling and, even when the detailed fits are used, it is difficult 
to distinguish between them. By combining the fits to the 
lines and the continuum around them, overall the best fit is 
for the OSP plus a 0.05 Gy r YSP. This result is entirely con- 
sistent with the findings of I Wills et aLl(|2004l ). although our 
results require a significantly lar ger contribution of the YSP 
(c.f. 10-30% in lWills et aLlbOO^I, and tow ards the older end 
of the range reported bv lAretxaea et al.l 12001!). 

Including a power-law component dramatically in- 
creases the number of models providing a good fit, 
marginally improves the Xred shifts the minimum to 
slightly older ages; OSP + YSP (0.1 Gyr, E(B-V) = 0.2, 
47% contribution) . Hence, for the nuclear region of Hydra A 
we find the best fitting model to be an old (12.5 Gyr) stellar 
population (45%) with a YSP (0.05 Gyr, E(B-V) = 0.4, 50% 
cont ribution). Th is is entirely consistent with the findings of 
both I Wills et all ([2004) (the difference in fiux percentage is 
likely to be due to differences in the seeing and extraction 
ap erture) and towards t h e older end of the range reported 
bv lAretxaga et all (|200ll '). IWills et al.l l|2004h also attempted 
to model the spectrum with younger age 'old' stellar popula- 
tions. They found the best fit comprised a 2 Gyr population 
(89%) with a 0.05 Gyr YSP (11%). We have also attempted 
to model the nuclear spectrum of Hydra A with 'younger' 
old populations and found this made no difference to our 
overall result. 



4.3.2 3C 236 

The radio source 3C 236 and its host galaxy are peculiar at 
both optical and radio wavelengths. The radio source exists 
both on the large scale (a classical double source extends 
to ~ 4 Mpc, making i t the largest known radio source, 
e.g. ISchilizzi et al.ll2001^ and yet the nuclear regions har- 
bour a second, compact (2 kpc) radio source suggesting it 
belongs to the cla ss of 'double-double' radio sources (e.g. 
iKaiser et al.l I2OOOI ) . On the basis of radio spectral ageing 
techniques and the detection of different star forming re- 
gions, two of which are relatively young (of order lO'^ yr) and 
two which are significantly older (of order 10^- 10^ yr, com- 



arab le to the age of the giant radio source*) . lO'Dea et al.l 
200 ll ) argue that this strange radio morphology is the re- 
sult of a period of fuel starvation of the AGN which lasted 
around 10^ years before new fuel was provided and the ra- 
dio source restarted. Note, the relative locations of the radio 
jets and the star forming regions preclude jet-induced star 
formation, but the timescales are consistent with a common 
origin of fuel such as the infall of gas to the host galaxy 
dO'Dea et al.ll200lh . 

At optical wavelengths the deep ground-based images 
of Smith & Heckman (1989) reveal the presence of shells 
and twisted isophotes that are suggestive of a past merger 
or interaction with a companion galaxy. On smaller scales, 
HST imaging shows the isophotes of the host galaxy to be 
flattened and inclined, with a broad asymmetric dust lane 
on scales of a few kpc which is itself misaligned w ith the 
inner disk (e.g. lMartel et al.ll 19991 : Ide Koff et al]|2000l). They 
also show a second, ^1 kpc scale dust feature ( de Koff et al.l 
2000.) . In the UV, three clear knots of emission are detected 
which trac e an arc coinciding with the SE edge of the larger 
dust lane (|Allen et al.ll2002l : lO'Dea et al.ll200lf ). 

The integrated nuclear spectrum of 3C 236 has a clear 
UV excess (Figure [1}. We find that the overall SED and the 
stellar absorption features can be satisfactorily reproduced 
(Xrod = 1-1) with an old (12.5 Gyr) stellar population (87%) 
and a red power-law component (a = 3.6, 12% contribution) 
and this model is shown in Figure!?] 

The power-law component required in this model is red. 
Following our 'decision tree' (Figure [2]), a red power-law 
component is unlikely to be a scattered AGN component 
and is more likely to be either direct AGN lig ht or a young 
reddened YSP component. The HST images o f Martel et al.l 
(| 19991 ): Ide Koff et all (|2000l ) and lO'Dea etal] l|200ll ) do not 
show a point source in the nucleus. Hence a large direct 
AGN light contribution is unlikely. Further, after subtract- 
ing the best-fit model we see no evidence for broad emission 
line components also expected if we are viewing the AGN 
directly. Hence, the UV excess is likely to be dominated by 
a very young YSP component which is heavily reddened. 

Figure [3] shows the x?od space for an OSP plus a YSP 
component. Two clear minima are observed for for YSPs 



with ages 0.1 Gyr (E(B-V) = 1.2, xi 



0.58) and 1.4 



Gyr (E(B-V) = 0.5, Xrcd = 0.48), although viable models 
(x?od < 1-0) to the SED are found across all ages. Detailed 
comparisons with the stellar absorption features rule out 
models with ages >1.0 Gyr because such models overpredict 
the depth of Call K feature, although it is difficult to dis- 
criminate further. An example of a two component model is 
shown in Figureg] 12.5 Gyr (84%) plus 0.05 Gyr YSP (15%, 
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E(B-V) — 1.4) with Xred = 0.78. Including a power-law com- 
ponent expands the number of viable fits. In summary, for 
the OSP plus YSP models, both with and without a power- 
law component, the YSP contributes typicaUy 10-30% of the 
flux in the normalising bin for YSP ages < 1.0 Gyr. Hence, 
our spectral fitting results are consisten t with both age esti - 
mates derived from HST photometry of lO'Dea et"aiT(|200ll ') 
- ~ 10^ yr and ^ 10*~^ yr - for different bright knots in the 
halo. 



4.3.3 3C 285 

The double-lobed FRII radio galaxy 3C 285 is identified 
wit h an elliptical g alaxy, the brightest in a group of galax- 
ies jSandagel I1972I ). A chaotic morphology consisting of a 
distorted S-shaped en velope aligned with another galaxy 
^40 arcsec to the NW (iHeckman et al.ll 19861 : iRoche fc Ealej 
I2OOOI ). tidal tails and fans, ir regular dust lanes and knots 
of emission (| Allen et al.l 120021 ) and a, B — V colour signifi - 
cantly bluer than a normal elliptical galaxy (|Sandagelll97^ ) 
all provide evidence that the galaxy has undergone a recent 
merger or interaction with a companion. The UV/optical 
continuum of 3C 28 5 has been previously modelled by 
lAretxaga et all (|200ll '). By fitting the 4000A break, they con- 
clude the blue light in the inner 2 arcsec is dominated by 
a sta rburst component with an age of 10-12 Myr. In addi- 
tion, Ivan Breugel fc Devi (|l993h have investigated the knot 
of UV emission located halfway between the nucleus and 
the eastern radio lobe and found the UBV colours and the 
4000 A break to be consistent with a 70 Myr starburst, that 
may have been triggered by the passing radio jet. 

The nuclear spectrum of 3C 285 is highly indicative of 
a major contribution to the continuum flux by light from 
young stars - the Balmer sequence is clearly observed in 
absorption and the galaxy displays a clear UV excess (Figure 
[T]). Indeed, SED modelling excluding a YSP contribution 
failed to provide an adequate fit to the overall shape of the 
continuum spectrum. 

Whilst it is necessary to have a YSP component to 
model the SED, good fits are obtained both with and with- 
out an additional power-law component. However, the emis- 
sion lines are relatively weak and HST imaging does not 
reveal a strong point sourc e; rather, the nucleus is obscured 
by o ne of the dust lanes (|Allen et all I2OO2I : iMadrid et all 
bOOej ). Models including a power-law can therefore be ex- 
cluded for two reasons: i) models requiring a very young 
YSP do not require a power-law component and ii) models 
with an older age YSP require a large power-law contribu- 
tion. We have therefore focussed on models which exclude a 
power-law component. The contour plot for the Xicd space 
for a 12.5 Gyr OSP and a YSP is shown in Figure [tj 

3C 285 is an ideal example to highlight the degeneracy 
in the solutions when relying on SED modelling alone - Fig- 
ure [3] shows two clear minima, one with a young (0.04-0.2 
Gyr), reddened (E(B-V) < 0.5) YSP contributing 30-50% of 
the flux in the normalising bin and another for a signiflcantly 
older (1.0-3.0 Gyr), unreddened YSP (E(B-V) < 0.2) con- 
tributing ~100% of the flux in the normalising bin. Detailed 
comparisons reveal that, whilst the fits to the majority of 
the diagnostic lines are indistinguishable, the viable models 
with an older YSP (> 1.0 Gyr) significantly overpredict the 
strength of the Call K line. Further, models with YSPs < 



0.1 Gyr do not fit either the absorption lines or the con- 
tinuum in the 3700-4300A region or the G band as well as 
other models. 

Hence, by combining both techniques, it is likely that 
the nuclear region of 3C 285 contains a significant YSP with 
age 0.1-0.5 Gyr and reddening 0.2 < E(B-V) <0.0, con- 
tributing 33%-39% of the continuum flux in the normalising 
bin. An example of a good flt is shown in Figure |4l 12.5 Gyr 
(72%) + 0.2 Gyr with E(B-V) =0.2 (33%); x?cd = 0.74. The 
age we estima te is at least an o r der o f magnitude larger than 
that found bv lAretxaga et all (120011 ) using the Balmer and 
4OOOA breaks and line indices (10-12 Mvr) . lAretxaga et al.l 
('2001') give no indication of the scale of the YSP contribu- 
tion to the flux. 



4.3.4 3C 321 

The powerful FRII radio source 3C 321 is identified with 
a p eculiar galaxy with two nuclei separated by ^3 arc- 
sec ( Spinrad et al.l Il985l : iHeckman et al.l 1 19861 : iBaum et al.l 
ll98S ) surrounded by an extended halo of line emission 
( McCarthy eTal] 1 19951 ). as well as continuum arcs and 
tails suggestive of a recent merger (IHeckman et al.l Il986l : 
ISmith fc HeckmanI Il989l : iRoche fc Ealesll2000l ). Similar to 
a significant number of powerful radio sources, 3C 321 
is blue in colour compa red to normal elliptical galaxies 
(jSmith fc Heckman|[l989l ). Figure [1] also highlights the spa- 
tially resolved UV excess. 

The SE nucleus is brighter and more closely coinci- 
dent with the radio c ore whilst the NW nucleus is clumpy 
(|Zirbel fc Baumlll998l ) but also associated with a peak in 
the radio emission which then extends as a lon g (35 arcsec) 
jet beyond this nucleus (|Roche fc EalesI I2OOOI ) . However, it 
is interesting to note that the narrow emission line spec- 
tra of the tw o nuclei have sim i lar line strengths and surface 
brightnesses (|Filippenkolll987l : iMcCarthv et al.lll995l ). 

3C 321 is a case in which the scattered quasar com- 
ponent makes a major contribu tion to the optical contin- 
uum, at least in the SE nucleus (iDraper et al.l|T99 ^). Spec- 
tropolarimetry observations of the SE nucl eus reveal a high 
degre e of polarisation in the near UV (|Tadhunter et all 
Il996l ). as well as broad quas ar-like permitted li nes in the 
polarised intensi t y spec trum ( Cohen et al.lll999|). H owever 



iTadhunte r et all (|l996l ) and 



Robinson et al 



X2OOO) have 



demonstrated that it is not possible to model the coritirmum 
SED solely in terms of an OSP plus a power-law or quasar 
template, but that a significant contribution from a YSP is 
required to fit the spectrum, particularly in the region of the 
Balmer break. A l thoug h th ey detected a YSP com ponent, 
ITadhunter et all (|l996l ) and (Robinson et all (|2000l ) consid- 
ered only unreddened YSP models, and were unable to pin 
down the age of the YSP. 

3C 321 is spatially well-resolved although, in contrast 
to lRobinson et al] (|2000l ). we have extracted only four aper- 
tures. These are shown in Figure [TJ which also highlights 
the significant UV exc ess observed acr o ss the entire galaxy. 
Despite the fact that iRobinson et al.l (I2OOOI ) have already 
presented work on these spectra, we have re-calibrated the 
data to improve the fiux calibration an d our analysis uses 
the la test high resolution models from iBruzual fc CharlotI 
(|2003l ) as well as consideration of the reddening of the YSP 
component. 
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3C 218 (Hydra A): 12.5 Gyr OSP + 0.05 Gyr YSP (E(B-V) = 0.4) 
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3C 285: 12.5 Gyr OSP + 0.2 Gyr YSP (E(B-V) = 0.2) 
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Figure 4. SEDs and detailed fits. Left: overall fit to SED where the black line traces the data, the red line is the overall model 
with one or more components (OSP, YSP, power-law) represented by the green lines. The wavelength bins used to evaluate the fit 
are identified by the purple boxes (not to scale) and the histogram gives the residuals in each bin. Right: detailed fits in the regions 
of the Balmer and Ca II lines. The solid line traces the data and the dotted line represents the model (red line in the left hand 
plots). 
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The SE nucleus. The spectrum of the SE nucleus of 
3C 321 shows evidence for the Balmer sequence in absorp- 
tion and so a YSP contribution is likely to be important. 
Indeed, previous studies of the continuum emision of 3C 321 
(see above) have shown that at least three distinct compo- 
nents are required to model the SED: an OSP, a YSP and a 
power-law component (the latter representing the scattered 
quasar indicated by the polarimetry measurements). Note 
that we cannot obtain viable fits with an OSP and power- 
law alone (x?cd> 2) in this aperture, consistent with the 
results of Tadhunter et al. (1996). Hence, we have focussed 
on the three component fits to the SED. 

When including a power-law component, as for many of 
the sources studied, acceptable fits to the SED were found 
for most ages of YSP - see Figure|31 Hence, the detailed com- 
parisons proved vital for distinguishing between the models. 
Using the Call K line (Figure|4]), we can confidently rule out 
models including a YSP component older than ~ 1.0 Gyr 

- all models with a YSP age > 1.0 Gyr significantly over 
estimate the depth of the Ca K line - this can not be over- 
come by smoothing of the models to match the resolution of 
the data. The youngest ages (0.04, 0.05 Gyr) are also con- 
sidered unlikely as they do not require a contribution from 
the power-law component (i.e. 0% contribution in the fits) 

- previous polarimetric observations combined with the de- 
tection of a broad component to Ha suggest a contribution 
of at least 10-20% at the position of our normalising bin 
(see above). Hence, the SE nucleus in 3C 321 can be well 
modelled by an OSP (15-55%) plus YSP (age: 0.1-1.0 Gyr, 
E(B-V): 0.1-0.3, 30-60%) and power law (10-20%) compo- 
nents with ~0 < a < ~ 0.6. Figure |4] shows an example of a 
model in the best fitting range comprising a 12.5 Gyr OSP 
(46%), a 0.4 Gyr YSP (E(B-V) = 0.2, 30%) and a power- law 
compon ent (a — 0.1, 22% ). These results arc consisten t with 
those of iTadhunter et al.l (|l996) and Robinson ct alj (|2000l ') 
in terms of the contribution o f the power-law compon ent and 
the YSP contribution. Note, iRobinson et all (|2000l ') do not 
find evidence for a YSP in the SE nucleus, only in the more 
extended regions. However, our SE nuclear aperture is larger 
than that used by Robinso n et al. 

I (Eooo) and therefore in- 
cludes light from more extended regions which do have a 
YSP component. 

The NW nucleus. As discussed above, whilst the SE 
nucleus is that identified with the radio source, the second 
nucleus ^^3 arcsec to the NW also shows signs of activity - 
it is close to an enhancement in the radio emission along the 
jet and has an emission line spectrum similar in strength 
and ionisation state to the SE nucleus. Further, as demon- 
strated by Figure[l] its observed UV excess is similar to that 
observed in the SE nucleus. 

As for the SE nucleus, it was not possible to obtain a 
good fit with an OSP, with or without a power-law {Xicd ~ 
2.0), and so a YSP component is likely to be important. 

It is possible to model the SED of the NW nucleus with 
an OSP plus a YSP without a power-law. Viable fits (x?cd ^ 
1.2) are found for all ages of YSP with the familiar 2 minima 
degeneracy. The 'younger' minimum is for ages ~0.1 Gyr 
(E(B-V) < 0.3) whilst the second minimum is for much 
larger ages (> 1.0 Gyr). Around the 'younger' minimum, 
the best fits to the lines and the continuum are for ages 
0.04-0.2 Gyr (12.5 Gyr: 70-75%, YSP: 25-30%) - ages > 0.3 
Gyr do not fit the continuum in the region 4000-4300A well. 



Solutions for all ages of YSP can provide good fits to the 
G-band and Mg lb features. 

Unlike in other apertures we have not been able to dis- 
tinguish between the two minima. At first glance, the mod- 
els including the older YSPs (1.0-1.4 Gyr) appear to signifi- 
cantly overpredict the strength of Call K. However, smooth- 
ing of the models to better match the resolution of the data 
show that for this aperture, the models including older age 
YSPs can not be confidently ruled out 

The strong emission lines detected in the NW nucleus 
suggest that i t also harbou rs a powerful AGN. Moreover, 
the results of [Draper et al.l (1993) show that it is signifi- 
cantly polarised. We have therefore attempted to model the 
SED of the NW nucleus with a three component model. In 
common with the two component OSP plus YSP model, we 
find viable fits for most ages of the YSP (0.05-2 Gyr) with a 
minimum consistent with the older ages in the two compo- 
nent models at ~ 1 Gyr. The contour plot of the Xrcd space is 
shown in Figure[3] Detailed fits show that the models includ- 
ing the oldest age YSPs (> 1.6 Gyr) can be ruled out due 
to the large over-estimate of the Ca K line and the models 
including the youngest age YSPs (<0.1 Gyr) are ruled out 
as they do not include a zero power law component. The re- 
maining models, however, are indistinguishable. Hence, the 
viable fits have YSP ages 0.1-1.4 Gyr with 0.1 < E(B-V) < 
0.3 contributing 25-75% of the SED flux. An example of 
a good fitting model is shown in Figure [4] 12.5 Gyr OSP 
(52%) plus 0.6 Gyr YSP (E(B-V) = 0.2; 30%) plus a power 
law (15%; a = -0.1) with xL = 0.45. 

Extended apertures Because the U V excess is spa- 
tially resolved, we have followed the work of lRobinson et all 
(j 2000l ) and extracted a further two apertures, one between 
the two nuclei and one to th e south east of the SE nucleus, 
marked on Figure [U Indeed, iRobinson et~al] l|2000l ) showed 
that, in the case of 3C 321, the YSP component makes a 
relatively larger contribution off-nucleus than in the regions 
of the two nuclei. 

Between the nuclei. It is not possible to model the re- 
gion between the nuclei with an OSP or OSP plus power-law 
component (x?cd > 2). This is not surprising as the spectrum 
shows clear evidence for a Balmer break and Balmer absorp- 
tion features, suggesting a strong YSP component. For OSP 
plus YSP models, viable fits (x?ed ^ 1-2) can be obtained for 
ages 0.05-2.0 Gyr with E(B-V) < 0.4, dependent on age - 
younger YSPs require more reddening. There is also a clear 
degeneracy with two minima at ages ~ 0.1-0.2 Gyr and ~ 
1.0-2.0 Gyr. However, examination of the overall fit to the 
SED shows that whilst the value of x?od is low, two regions 
are not modeUed well - ~5000-5500A and >7000A. 

As the two component models fail to successfully fit 
these two regions, we have attempted to model the SED 
including a power law. Whilst for most off-nuclear apertures 
we have assumed that a power law component is unlikely, 
for 3C 321, based on the previous observational results, this 
is not necessarily the case. 

Figure [3] shows the x?od space for the three component 
models. When including a power law, the Xicd reduces sig- 
nificantly and a broad minimum is found for ages 0.04-2.0 
Gyr with reddening E(B-V) < 0.5. Combining the large scale 
and detailed fits, the youngest ages can be ruled out on the 
basis of the large scale (5000-5500A and >7000A regions) 
fits (0.04-0.1 Gyr) and the fits to the Call K line (< 0.3 
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Gyr). Older age YSPs (> 1.0 Gyr) can be ruled out as the 
Call K line is significantly over-predicted. 

Hence, the best fitting models for the region between 
the nuclei comprise a 12.5 Gyr OSP (30-50%) plus a 0.4- 
1.0 Gyr YSP (E(B-V) = 0.0-0.1, 40-60%) and a power law 
component (5-8% with a ~ -1.6 - -1.7). An example of a 
model in this range is shown in Figure [l] 12.5 Gyr plus 0.4 
Gyr YSP (E(B-V) = 0.1, 39%) plus a power law (8%, a = 
-1.6) with xL = 0.42. 

Extended aperture. As for the aperture between the 
nuclei, it is not possible to model the SED with an OSP or 
OSP plus power- law alone {x^ad > 6). ModeUing the SED 
using two and three components, each incorporating a YSP, 
shows that, whether or not a power-law component is in- 
cluded, the best fits to the overall SED consistently required 
a ~0.7 - 2.0 Gyr YSP with zero reddening. Figure |3] shows 
the Xred space for both two and three component models. 
This is reassuring and shows our results are robust against 
the fine details of the modelling. Detailed comparisons of 
the absorption lines show that, when a power-law compo- 
nent is excluded, the Call K line is always significantly over- 
predicted. Including a power-law component has the effect 
of diluting the absorption features. Indeed, this produces a 
significantly better fit to the Call K line for a YSP age of 
0.7-1.0 Gyr, whilst older and younger ages give worse fits 
to either the Call K line and/or the continuum surrounding 
the absorption features. 

The modelling suggests that three components are re- 
quired to model the spectrum. A significant YSP component 
(54-72% of the flux in the normalising bin) is essential to pro- 
vide a good fit to the SED and detailed comparisons reveal 
that a small (5-6% of the fiux in the normalising bin) power- 
law component improves the fit to the absorption features. 
An example of a good fitting model, 12.5 Gyr (21%) plus 
1.0 Gyr YSP (E(B-V)= 0.0; 72%) plus power-law (5%; a = 
-1.2) with Xred = 0.7 is shown in Figure 3] This aperture 
contains the largest contribution from a YSP (>50% in the 
normalising bin). 

Summary for 3C 321. The good spatial resolution of 
the continuum structure in 3C 321 allows us to investigate 
the spatially extended UV excess in four distinct regions in- 
cluding both nuclei, the region between the nuclei and an 
extended aperture to the SE of the active nucleus. Despite 
the large spatial scale and the different environments anal- 
ysed, the results are generally consistent. All apertures re- 
quire a significant contribution (~30-70%) from a ~ 0.5-1.0 
Gyr YSP with little reddening (E(B-V)< 0.3) although we 
can not rule out significantly younger ages in the two nuclei. 
In addition, all apertures require a relatively small contribu- 
tion (few-20%) from a power law in all apertures, the largest 
contribution in the active nucleus. 



4.3.5 3C 381 

The ultra-steep spectrum radio source 3C 381 is a clas- 
sical double r adio source extending more than 70 arcsec 
along PA 90 (|Rilev fc PoolevI Il975l i. Optically it is iden- 
tified with a round high surface brightness galaxy with 
a compact nucleus. A lower surface brightness compan- 
ion lies 5.5 arcsec away along PA 65 and the detection 
of tidal distortions and tails extending N and S con- 
firms 3C 381 is currently interacting with this compan- 
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Figure 5. Detection of a broad component to Ha in the nuclear 
aperture of 3C 381 with FWZI ~ SOOA or ~ 14,000 km s"!. The 
lines plotted are: nebular continuum subtracted nuclear spectrum 
(bold solid Une; offset by -0.4 X 10~^^ erg s~l cm~^ best 
fitting (OSP+power law; offset by -0.4 X 10"-*^^ erg cm~^ 
A~^) SED model (solid line) and the continuum (OSP-|-power 
law) subtracted nuclear spectrum (dotted line). Note, the spec- 
trum has not been corrected for atmospheric absorption features 
and the A band (at ~ 7600 A) Is coincident with the Ha/[N II] 
blend In the observed frame. 



ion jd e Koff ct al' 'l996"; 'Roche fc Eaie3 |2000| ). Ahhough 
ICrand i & Ostcrbrock (1978) reported the detection of weak 
broad wings to Ha along wit h high ionization narrow lines , 
iDennett-Thorpe et al.l (I2OOOI ') fc lEracleous fc HalpenI (120041 ) 
failed to find ev idence for broad Ha in t heir h igher S /N spec- 
tra. Moreover, iDennett-Thorpe et al.1 l|2000l ) argue against 
the classification of this object as a BLRG because the 
radio core and jet are relatively weak, and the near-IR 
colours are consistent with NLRGs r ather than BLRGs (c.f. 
iLillv fc Longaiilll982D . iRudv et al.l (|l983D measure a rela- 
tively low polarisation for 3C 381 at optical wavelengths: 
1.8±0.6%, consistent with the observations of other BLRGs 
in their sample. On the basis of Figure[T] 3C 381 has a clear 
UV excess. Given that it is only marginally resolved (Figure 
[1}, we have extracted only a nuclear aperture for this source. 

We find that an adequate fit to the SED and the discrete 
features can be obtained using a two component (OSP plus 
power-law) model (x?cd = 0.90) where 80% and 20% of the 
flux originates from the OSP and power-law respectively. 
This model is shown in Figure 3] In this model, the power- 
law component is red (a = 1.65) which suggests that, if 
the light originates from the AGN, it is direct rather than 
scattered, consistent with the low optical polarisation. 

In Figures |4] and [5] there is clear evidence for a broad 
component to Ha in the form of an excess in the continuum 
region around Ha relative to our best fitting model. This 
excess has a FWZI of ~300A or ~14,000 km s"\ However, 
due to the strong atmospheric A band (~ 7600A) coinci- 
dent with the Hq/[N H] blend, it is not possible to model 
this component i n det ail. We further note that, although 
IChiaberge eT^ (|2002l ) present only an upper limit on the 
flux of a compact nucleus in this source based on HST im- 
ages, their upper limiting flux is consistent with the flux in 
the power-law component required by our continuum mod- 
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3C 321 (SE nuc): 12.5 Gyr OSP + 0.4 Gyr YSP (E(B-V) = 0.2) + power law 
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3C 321 (NW nuc): 12.5 Gyr OSP + 0.6 Gyr YSP (E(B-V) = 0.2) + power law 
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3C 321 (mid, scaled): 12.5 Gyr OSP + 0.4 Gyr YSP (E(B-V) = 0.1) + power law 
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3C 321 (ext): 12.5 Gyr OSP + 1.0 Gyr YSP (E(B-V) = 0.0) + power law 
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Figure 4. SEDs and detailed fits continued 



elling. Hence, our observations are consistent with a BLRG 
classification for 3C 381. 

Despite tfie possible evidence for a broad line nucleus 
in this object, the evidence is not conclusive. Therefore, 
we have modelled the SED including a YSP, both with 
and without a power-law component. Figure [3] shows the 
X?od space for the OSP plus YSP models. For the two com- 



ponent models, we find one mininmum (Xrod = 0.85) for an 
OSP plus a 0.02 Gyr YSP reddned by E(B-V) = 1.0 con- 
tributing 82% and 17% of the flux in the normalising bin 
respectively. This model also produces a viable fit to the 
stellar absorption lines. 
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4.3.6 3C 433 

The FR I/FR II radio galaxy 3C 433 is peculiar at both ra- 
dio and optical wavelengths. The radio source extends over 
~ 100 kpc and is highly asymmetric in both the strength and 
shap e of the radio emission, overall resembling an 'x' shape 
(e.g. Ivan Breugel et al.lll983l : iBlack et al.1ll99^ ). Such com- 
plex radio morphology, often associated with lower power 
radio sources, is une xpected due to the relati vely high power 

of the radio source (jvan Breugel et al]ll983l ). 

^Optically, 3C 433 is identified with a NLRG l|Koskil 

Il978l), with a bright compan i on galaxy ~10 arc sec to the 
NE Ivan Breugel et al.1 [l983 : iBaum et all 119881 ) . and the 
nuclei of both galaxies are surrou nded by a common en- 
velope (|Matthews fc Morgan|[l963 ). A fainter galaxy is ob- 
served ~25 arcsec to the SW. The large-scale optical struc- 
ture is complex, with emission line fil aments clearly de- 
tected in grou nd-based (Ha) images (e.g. iBaum et al.]|l988l : 
iMcCarthv et a l. 1995), and patchy dust features, some trail- 
ing in the direction of t he radio jet, also de tected in higher 
resolution HST images l|de Kofi et al.ll200Gl ). 

Modell i ng of the spectrum of the SW nucleus by 
IWills et al.l ([2002!) suggests that the large UV excess ob- 
served is attributable to a significant (~ 20-40% contribution 
to the flux at 4780A), highly reddened (E(B-V) ~ 0.5-0.7) 
YSP (0.05-0.1 Gyr). Note that, on the basis of the lack of a 
broad component to Ha and the low UV polarisation (< 8%; 
I Wills et al.|[2003 ). reddened YSP models are favoured over 
scattered AGN models for this source. Moreover, while this 
object shows unu sually red colours at near-IR wavelengths 
iLillv et al.ll 19851 ). and a relatively high surface brightness in 
the complex nuclear structure revealed by the optical HST 
images, the nucleus is resolved and there is no evidence for 
a point source at optical wavelengths (Chiaberge, Capetti & 
Celotti 1999). 

In our new WHT observations the slit position was cho- 
sen to cover both nuclei. We find that the UV excess covers 
the entire measureable extent of the double nucleus system 
(see Figure[l|. Hence, we have extracted three apertures: one 
for the SW nucleus coincident with the radio source, one for 
the NE nucleus, and one sampling the region between the 
nuclei. All apertures are shown on Figure [1] 

The SW nucleus. The SW nucleus of the 3C 433 sys- 
tem, coincident with the main radio source, is likely to have 
a YSP component (see Wills et al. 2002). However, it is pos- 
sible to obtain a marginally acceptable fit (Xrcd ~ 1-40) to 
the SED and absorption features using an OSP (74%) plus a 
power-law component (24%; a = 1.0) only, cons istent with 
the two component fitting by I Wills et al.l ( 2002 ). Since the 



UV polarisation is low (<8%; Wills et all 



20021 ). the large 



power-law is unlikely to be due to a scattered AGN com- 
ponent. Although the fitted power-law is red, which might 
be consistent with a direct AGN component, given the fail- 
ure to detect a point source in broad-band HST images 
l|Chiaberge et al.]|l999l ) and broad wings to the Ha in our 
spectra, an interpretation in terms of a highly reddened YSP 
component is more plausible. 

Figure [3] shows the plot of the x?od space for models 
that o nly include YSP and OSP components. IWiUs et al.l 
l|2002D found a two component OSP plus YSP model gave 
as good or better fits than an OSP plus a power-law and the 
results presented here are consistent with this - there is a 



clear mininmum in the Xrcd space at YSP age ~ 0.05 Gyr 
with reddening E(B-V) ~ 0.7 and Xrcd ~ 0.6 contributing 
~ 20-30% of the flux in the normalising bin. Models with 
YSP ages greater than 0.2 Gyr are ruled out by our SED 
modelling. An example of a two component model is shown 
in Figure m with: 12.5 Gyr OSP (67%) plus 0.05 Gyr YSP 
(E(B-V) = 0.7, 32%) with x?cd = 0.6. It is clear that this 
model also provides an acceptable fit to the stellar absorp- 
tion features. 

The NE nucleus. The SED of the NE nucleus is 
markedly different to that of the active SW nucleus - the 
4000A break is more prominent, as are the stellar absorp- 
tion features (e.g. Call H-l-K) , whilst there is little evidence 
for line emission. However, it is possible to model the SED 
with an OSP, both with and without a power law (x?od = 
0.99 and 1.1 respectively). Figure g] shows the OSP (100%) 
model. 

Despite being able to model the SED adequately with- 
out a YSP component, because there is evidence for a signif- 
icant UV excess (see Figure [1}, we have attempted to model 
the SED including a YSP to see how the fit will change. 
In the two component OSP plus YSP models, the SED fit 
improves significantly when including a YSP and good fits 
iXrod < 1) S're obtained for all ages of YSP (see Figure 
The inclusion of the YSP component also improves the fit 
to the absorption lines, although the younger (< 0.03) and 
older (> 1.6) ages tend to over predict some of the absorp- 
tion features. An example of the 12.5 Gyr OSP (89%) plus 
a 0.1 Gyr YSP (E(B-V) = 0.8, 10%) with x?ed = 0.57 is 
shown in Figure U) Despite the fact that the models are vir- 
tually indistinguishable, all consistently have OSP and YSP 
contributions of 95-65% and 5-30% respectively with x?cd ~ 
0.45-0.9. 

Between the nuclei. As highlighted in Figure [T] the 
UV excess is spatially resolved and extends across the entire 
system, with a similar _D4000 ratio. We have therefore also 
investigated an aperture located between the two nuclei. The 
location of this aperture is shown on Figure [T] 

It is not possible to model this aperture with an OSP 
alone, with or without a power-law component. Moreover, 
an AGN component is not expected in this aperture as it 
is located far from the active nucleus and does not contain 
any line emission - the UV excess is therefore likely to be 
due to a population of young stars. 

Figure [3] shows the x?cd space for two component OSP 
plus YSP models. Two minima are clearly seen at ~ 0.1- 
0.2 Gyr (E(B-V) ~ 0.5) and ~ 1-2 Gyr (E(B-V) ~ 0.2) - 
note that viable fits are also found between the minima. 
As for other apertures displaying such a degeneracy in the 
fits, we can rule out older age YSPs on the basis of the 
strength of Call K - models including YSPs older than ~ 
1.0 Gyr significantly over-predict the strength of Call K. For 
all models in this range, the OSP and YSP contributions are 
typically 80-95% and 10-25% respectively with E(B-V) ~ 0- 
0.7. Figure |4] shows an example of a good fitting model: 12.5 
Gyr OSP (91%) plus 0.1 Gyr YSP (E(B-V) = 0.5, 15%) with 
xL = 0.95. 



4.3.7 PKS 0023-26 

The powerful (log Psghz = 27.33 W Hz"^) Compact 
Steep Spectrum (CSS) radio source PKS 0023-26 is iden- 
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Figure 4. SEDs and detailed fits continued 



tified with an nit, = 1 9-5 galaxy JPrestage fc PeacocklflQS^ : 
IWall fc Peacocklll985l ). The double lobed radio source has 
a maxi mum size of ^ 680 m as (~ 4.25 kpc), aligned along 
PA -34 (|Tzioumis et al.ll2002l ). Spectroscopically it is classi- 
fied as an NLRG and its emission line spectrum has inter- 
mediate ionisation w ith relatively str ong [O n]AA3727 and 
[N l]A5f99 emission (jTadhunter et al.iif99a : .Morganti et all 



If997f ). The optical continuum displays a large UV excess 
with low UV polarisation {Puv < 2.0%) and no sign of 
broad permitted lin es. This source is therefore a possible 
starburst candidate ('Ridhunter et al.|[l993 : iMorganti et al.l 
ll997l : lTadhunter et al.ll200^ 

Figure [T] shows that the UV excess is extended beyond 
the immediate nuclear regions (indicated in figure). We were 
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unable to obtain adequate fits to the spectrum of the nuclear 
aperture of PKS 0023-26 with an OSP alone or an OSP plus 
power-law component. Indeed, the spectrum shows clear ev- 
idence for the high order Balmer lines in absorption (see 
Figure |4l 

For models including an OSP and a YSP we find a clear 
minimum in the xLd space for a young YSP age (~ 0.03-0.05 
Gyr) with significant reddening (E(B-V) ~ 0.8) and a ~ 50- 
60% contribution from the YSP component (see Figure |3]). 
Detailed comparisons for this aperture are difficult as many 
of the absorption features are infilled by strong emission 
lines. However, the Call K line is not filled and on the basis 
of the fit to Call K and the surrounding continuum, the 
best fitting model is for an OSP plus a young (0.03 Gyr) 
YSP with significant reddening (E(B-V) = 0.9); this model 
is shown in Figure H 12.5 Gyr (44%) plus YSP (0.03 Gyr, 
E(B-V) = 0.9, 54%) with x?od = l-l- Models including older 
age YSPs fail to reproduce the continuum in this region and, 
as the age increases, the Call K line is over-estimated. 

As for many of the other objects/apertures, including 
a power-law component significantly increases the range of 
viable YSP fits. A single broad minimum is observed cover- 
ing a large range of YSP ages (0.05-2.0 Gyr) and reddenings 
(E(B-V) < 1.2). YSP models within this minimum are in- 
distinguishable due to the power law contribution (10-60%). 
However, previous studies of PKS 0023-26 have found that, 
whilst the UV excess is large, there is little evidence for 
a significant AGN component - no broad permitted lines 
are o bserved (making a sign ificant direct component un- 
likely, iTadhunter et al.1 12OO2I and this paper) and the UV 
polarisation is low, ruh ng out a large scattered component 
ITadhunter et al.ll2002D . 

It is encouraging to note that all of the better fits includ- 
ing a power-law component require at least 20% of the flux 
to originate from the YSP, often as high as 50%, although 
the dominance of the power law component for some models 
make distinguishing between the ages in the three compo- 
nent models difficult. 



4.3.8 PKS 0039-44 

PKS 0039-44 is one of the higher redshift sources in our 
sample {z — 0.346) and has been previousl y studied as part 
of the 2 Jy survey by Tadhunter et al. (e.g. IXadhunter et al.] 
I1993I . I2OO2I '). This FR-II double-lobed radio source is identi- 
fied with a gal axy displaying a r ich, high ionisation emission 
line spectrum (|Tadhunter et al1 [l993l and is classified as an 
NLRG. 

PKS 0039-44 shows evidence for a s i gnifica nt UV ex- 
cess (see Figure [T]) and ITadhunter et al.] l|2002l ) have pre- 
viously investigated the optical/UV continuum using lower 
resolution data and spectral synthesis models. Whilst their 
observations show PKS 0039-44 to be significantly polarised 
{Puv = 4.8±1.2%), sug gesting a scattered AGN component 
is likely to be significant. lTadhunter et al.l l|2002l ) believe this 
component is not necessarily the dominant component in 
the UV. This is because the intrinsic polarisation of the 
continuum component causing the UV excess is relatively 
low {Puv < 10%), and they detected no broad emission 
lines. Tadhunter et al. therefore believe that PKS 0039-44 is 
a good candiate for a significant contribution from a YSP. 



ITadhunter et al.] (|2002l ) also report a 19% contribution of 
the nebular continuum in the UV. 

As found bv lTadhunter et al.1 l|2002l ). we can adequately 
model the SED using an OSP and a power-law component 
(Xrcd — 0.4) contributing 57 and 40% to the fiux in the nor- 
malising bin respectively. As well as reproducing the SED, 
this fit also provides an adequate fit to the stellar absorption 
features, and the slope of the required power-law (a ~ —0.5) 
is consistent with a scattered AGN component. Further ev- 
idence for a scattered AGN component is provided by the 
clear detection of a broad Mgll feature (FWHM ~ 9000 km 
s~^: see Figure [6)|. 

Although the detection of significant UV polarisation 
and a broad Mgll feature suggest a major contribution from 
a scattered AGN component, this does not by itself rule out 
a significant contribution from a YSP component. Hence, we 
have modelled the SED including both YSP and power-law 
components. Unsurprisingly, given the contribution of the 
scattered AGN, the models do not provide strong constraints 
on the age of the YSP: all ages of YSP provide a viable fit to 
the SED (Xrcd IS 0.4) and the detailed fits to the absorption 
lines are indistinguishable. However, by comparing the per- 
centage contributions of the components, some models can 
be ruled out if it is assumed that both a YSP and a power- 
law are required - many of the viable fits do not require a 
YSP component (0% model contribution). Ruling these out, 
all of the results are broadly consistent with component con- 
tributions of ~ 20-50%, ~ few-20% and 25-40% for the OSP, 
YSP and power-law (-0.7 < a < -0.2) respectively. However, 
the this only rules out certain regions of the reddening space 
(E(B-V) < 1.0 with larger reddenings required for younger 
YSPs) and viable models are still obtained for all ages of 
YSP. Overall, while the optical continuum of PKS 0039-44 
clearly has a large contribution from a scattered AGN com- 
ponent, we cannot rule out a significant contribution from 
a YSP conponent. Indeed, such a component would help to 
explain the fact that the UV polarisation is lower than ex- 
pected if the UV excess is solely due to scattered AGN light 
ITadhunter et al.ll2002l '). 

PKS 0039-44 is one of the higher redshift objects in our 
sample {z = 0.346). Therefore we have further modelled the 
SED using different OSPs in addition to the 'standard' 12.5 
Gyr population used for the majority of the sources in the 
sample. Based on the redshift, a Hubble constant of 75 km 
s~^ Mpc~^ gives ages of the universe at the source redshift 
between 10.5 and 13 Gyr for an open and ffat universe re- 
spectively. We have therefore also modelled the data with 
an OSP of age 10 Gyr. This gives similar results (50-55% 
old, 40-45% power-law with 0.4 < a < 0.5) to models with 
a 12.5 Gyr OSP. Therefore we conclude that, given the large 
contribution of the power-law and/or YSP components, our 
modelling does not provide strong constraints on the age of 
the OSP component. 

Because the SED of PKS 0039-44 can be adequately 
modelled by an OSP and a power-law component and there 
is the hint of a possible broad component to Mg II, we do 
not show a contour plot for the YSP models. As an exam- 
ple of the fits obtained, we show the model comprising an 
OSP and a power-law component with the best chi-squared, 
comprising a lOGyr OSP (57%) and a power-law component 
(41%; a = -0.5) in Figure H 
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Figure 4. SEDs and detailed fits continued 



4.3.9 PKS 0409-75 

At 2 = 0.693, PKS 0409-75 is both the most distant ra- 
dio source in our sample and the most powerful at radio 
wavelengths. Indeed, it is one of the most powerful inter- 
mediate redshift radio sources in the southern hemis phere 
(L45MHZ-4.8GHZ = lO^^'^^ erg s'^ lAlvarez et al.lll993l ). The 
FR II radio source is relatively small (LAS = 9 arcsec or 85 



kpc) and is domin ated by two lobe s of emission extending 
along PA130 (Mo rganti et all 1 19991 ). PKS 0409-75 is opti- 
cally identified with an my = 21.6 galaxy which is possibly 
resolved into two com ponents and displays extended emis- 
sion in [O II]AA3727 |di Serego Alighieri et all 1 19941 ). Low 
resolution optical spectra of PKS 0409-75 reveal a signifi- 
cant UV excess compared with passively evolving elliptical 
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galaxies at similar redshifts (|Tadliunter et al1l2002l ). This is 
also highlighted in our new measurements shown in Figure [1] 
The narrow emission line luminosity and ionisation state are 
relatively low for the radio power - [O III] is barely detected 
whilst [Q II] is re l atively strong in the spec t rum p resented by 
iTadhunter et all l|l993l i. Ffadhunter et all (|2002r ) attempted 
to model the optical/UV continuum but were unable to ob- 
tain a result due to poor sky subtraction. They do report a 
measurement of the 4000A break which suggests 95% of the 
light in in the 3750-3950A region is not emitted by the OSP. 
Further, due to the low UV polarisation {Puv < 6.3%) and 
non-detection of broad permitted lines in the UV spectrum 
(e.g. Mg II), Tadhunter et al. propose PKS 0409-75 to be a 
strong starburst candidate. 

By using the VLT and paying careful attention to the 
calibrations and sky subtraction, we have obtained new 
deep, higher resolution optical spectra of PKS 0409-75 which 
are of sufficient quality to model the SED. Using the simplest 
models, we find the SED shape can be adequately modelled 
using an OSP and a power-law component. As for PKS 0039- 
44, we have also considered OSP ages younger than 12.5 Gyr 
and find the best fitting model is for a 7 Gyr OSP plus a 
power-law component (x?ed = 1-2), each contributing 60 and 
45% of the flux in the normalising bin respectively. This flt 
is shown in Figure |4] 

However, whilst an adequate fit can be gained without 
a YSP component, no significant p olarisation has been de - 
tected for this source in the UV jTadhunter et all l2002h . 
the required power-law spectral index is relatively red 
(a ~0.86), no broad permitted lines have been detected, 
and the emission line luminosity is relatively low. Together, 
these features make a major contribution from a scattered 
AGN component unlikely, although without high resolution 
HST imaging we cannot absolutely rule out a contribution 
from a direct AGN component. 

Given the lack of clear evidence for an AGN contribu- 
tion to the UV excess, we have also attempted to model the 
SED of PKS 0409-75 including a YSP component. 

For models with an OSP (7 Gyr) and a YSP only, viable 
fits can be obtained with a young (< 0.04 Gyr), reddened 
(E(B-V) ~ 0.8) YSP contributing ~ 25-45% of the flux in 
the normalising bin. Figure |3] shows the Xrcd space for mod- 
els including an OSP and a YSP and Figure |4] shows an 
example of a good fitting two component model: 7 Gyr OSP 
(73%) plus 0.02 Gyr YSP (E(B-V) = 0.7; 32%) with xL = 
0.75). Inspection of the detailed fits shows that, due to the 
noise in the data, all of the models close to the Xied min- 
imum provide good fits to the absorption features and are 
indistinguishable. 

4.3.10 PKS 1932-46 

The double-lobed FR II radio source PKS 1932-46 is op- 
tically identified with an early-type galaxy at z = 0.231 
located b etween the two radio lobes but closest to the east- 
ern lobe (|di Serego Alighieri et all Il994h . The galaxy has 
a double structure in [O III] emission with the two com- 
ponents separated by 5 arcsec (25 kpc). The SW com- 
ponent is coincident with the optical continuum emission. 
The main galaxy structure is surrounded by a giant (up to 
~ 160 kpc) emission line nebula, making this one of the 
largest emission line haloes ever detected around an ac- 



tive g alaxy at any redshift lIViUar-Martm et al.|[l99i . 120051 ). 
Note, IVillar-Martm et all (120051 ) used the same data set 
presented in this paper. The halo is complex in structure 
with clumps and filaments of emission and evidence for stel- 
lar photoionization of the emission line regions on scales of 
more than 100 kpc in a direction close to perpendicular to 
the radio axis (|Vmar-Martm et allbOOSl). The nuclear s pec- 
trum shows a clear UV excess ( Tadhunter et all |2002| : see 
also Figure [TJ, a nd a rich emission line spectrum with mod- 
erate ionisation ([Tadhunter et al1ll993l : IVillar-Martfn et al.l 
l2005h . The detection of weak bro ad permitted lines (FyVIIM 
~ 2400 ki n s~M was reported bvlVillar-Martm et all (|l998l ) 
although IVillar-Martfn et al.l l|2005h retracted this detec- 
tion after analyses of higher quality data showed Hq could 
be modelled by the same kinematic c omponents as [O 
III]AA4959,5007. ITadhunter erall tOOj ) have modelled a 
lower resolution spectrum of the nuclear regions and sug- 
gest that, in addition to an OSP and power-law component, 
a contribution from a YSP is also likely. PKS 1 932-46 is 
not s i gnificantly polarised in the UV {Puv < 1.6%: iDlcksoni 
1 19971 : iTadhunter eral]|2002l ) . 

We find that it is possible to model the nuclear SED 
of PKS 1932-46 adequately with an OSP and a flat(a = 
0.23) power-law component (x?od = 0.82) with 59% and 41% 
of the fiux in the normalising bin originating in the two 
components respectively. The best fitting model is shown in 
Figure |4l 

Given that the UV polarisation is low l|Tadhunter et all 

and the modelled power-law component is red, this 
suggests that any light originating from the AGN will be 
direct. The best (ground-base d) images of PKS 1 932-46 
known to us are presented in IVillar-Martm et al.l (|l998l ) 
and are of insufficient resolution to determine whether a 
strong point source is detected. However, further evidence 
for or against a significant AGN contribution can be gained 
from the detection/non-detection of broad permitted lines. 
As discussed above, previous searches for broad lines in 
the nuclear spectrum give contradictory results. After sub- 
tracting off the best fitting model (12.5 Gyr OSP plus 
power- law), we find evidence for a broad component to 
Ha (see Figure [7} although this component is significantly 
broader (FWHM ~ 5100 km s~^) tha n the original detec- 
tion (FWHM ~ 2400 km s^^). Whilst IViUar-Martfn et al] 
( 2005) used the same spectra as presented here and searched 
for such a broad component, due to its broadness their non- 
detection is not surprising for a number of reasons. First, 
although the fiux in the component is relatively large, the 
large FWHM reduces the peak flux to the point that it is 
only clear when a large spectral ra n ge is studied. Second, 
analysis in the IVillar-Martm et al.l (|2005l ) paper focussed 
on a region around Ha only 250A wide in the observed 
frame whilst the component detected here has FWHM ~ 
140A and a total width of ~_370A in the observed frame. Fi- 
nally, IViflar-Martm et all (|2005h did not model and subtract 
the continuum - a step which has been shown to be vital 
for studying the broader emissi on line compone nts in radio 
galaxies (e.g. in PKS 1345-H2: iHolt et al.ll20o3 ). We there- 
fore support the classification of PKS 1932-46 as a BLRG. 

Given the strong evidence for broad lines in the nu- 
clear spectrum of PKS 1932-46 it is likely that a direct AGN 
component makes a major contribution to the optical/UV 
continuum. However, there is strong evidence from the HII- 
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Figure 6. Detection of a broad component to Mg II in the nuclear aperture of PKS 0039-44. a) Plot showing a zoomed in region 
of the nuclear spectrum of PKS 0039-44 highlighting the underlying broad component to Mg II. The lines plotted are: nebular 
continuum subtracted nuclear spectrum (bold solid line), best fitting (OSP+YSP-|-power law) SED model (solid line) and the 
continuum (OSP+YSP-|-power law) subtracted nuclear spectrum (dotted line). In plot b) we model the emission line components 
of the Mg II blend. The 2 component model comprises an intermediate (~ 1100 km s~^) component (dotted lines) and a broad (~ 
9300 km s~^) component (dashed lines). 
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Figure 7. Detection of a broad component to the emission lines in the nuclear aperture of PKS 1932-46. a) Plot showing a zoomed 
in region of the nuclear spectrum of PKS 1932-46 highlighting the underlying broad component to Ho. The lines plotted are: nebular 
continuum subtracted nuclear spectrum (bold solid line), best fitting (OSP-|-power-law) SED model (solid line) and the continuum 
(OSP+power-law) subtracted nuclear spectrum (dotted line). In plot b) we model the emission line components of the Ha/[N II] 
blend. The 7 component model comprises a narrow (180 km s~^) component (dotted lines) and a broad (1550 km s~^) component 
(dashed lines) to both Ha and the [N II] AA6548,6583 doublet, and a further very broad (5130 km s~^) component (dot-dashed 
line) to Ha. 



region like emission line ratios for ongoing star formation in 
the extended halo of this object (|Villar-Martm et al.ll2005l ) 
Moreover, when fitting a two comp onent model (OSP plus 
power- law) , Tadhunter et al.l (|2002t ) remained confident that 
a YSP component was required after the subtraction of their 
model left significant residuals. Indeed, detailed inspection 
of the SED in the region ~3700-4100A reveals a narrow ab- 
sorption hne at ~ 3830 A- most likely H9 A3835. Whilst the 
two component (OSP plus power-law) model reproduces the 
majority of the SED well, it fails to provide a good fit to 



this feature. Hence we have also attempted to model the 
SED using a three component fit (as the evidence suggests 
a significant AGN component is required) to include a YSP; 
the plot of the Xrcd space is shown in Figure |3l Unsurpris- 
ingly, the three component models produce viable fits for 
all ages of YSP with a large minimum in the range 0.05-3.0 
Gyr with E(B-V) <0.8 and it is virtually impossible to dis- 
tinguish between these models. However, it is clear that a 
YSP component is required to model the H9 line. The mod- 
els have varying contributions from the three components: 
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Figure 8. Identification of the features in the nuclear aperture 
of PKS 21 35-20 in the UV in the range 2200-3400A. Line IDs are 
taken from lTran et al.l lll998l ). The nebular continuum subtracted 
spectrum of the nuclear aperture of PKS 2135-20 is plotted along 
with a scaled version of the SDSS quasar template taken from 
IVanden Berk et al and line identifications for the known 

emission lines in this region. From this, it seems likely that the 
broad feature at ~ 2600-2700A is the result of several differ- 
ent emission line components, namely a broad component to Mg 
IIA2798 plus emission from He IIA2734 and [Mg VII]A2632. The 
dip observed in the spectrum of PKS 2135-20 at ~ 2600A, which 
is not observed in the quasar template spectrum, is possibly at 
least partially caused by the combination of the quasar with the 
light from the stellar populations — there is a small 'break' in all of 
the stellar templates at ~ 2600A present which appears stronger 
for older age populations. 



- 0-60% OSP, ~ 10-80% YSP and ~ 5-40% power-law with 
a > -0.6. Including the YSP improves the fit to H9, with the 
marginally better fits for the younger age YSPs suggesting 
that a YSP component may be important on the 10-20% 
level in the normalising bin. In addition to the OSP plus 
power-law fit, we also show an example of a three compo- 
nent model in FigureH 12.5 Gyr (54%) plus a 0.05 Gyr YSP 
(19%; E(B-V) = 0.4) and a power-law component (28%; a 



= 0.9). This model has a Xrcd of 0.44. 



4.3.11 PKS 2135-20 

The powerful (log Psghz = 27.47 W Hz"^) Compact Steep 
Spectrum (CSS) radio source PKS 2135-20 is identified with 
an mv = 19.4 galaxy (IWall fc Peacocklll985l ) at redshift z 
= 0.635 (|di Serego Alighieri et all 1994h and has the second 
highest redshift in our sample. At radio wavelengths, the 
MERLIN 5 GHz image shows a double-lobed, highly asym- 
metric structure cover ing ~250 mas (~3 kpc) along PA 52 
jTzioumis et all 120021 '). All of the radio fiux is detected in 
the MERLIN image suggesting there are no extended com- 
ponents. 

Broad-band imaging of PKS 2135-20 by |Sh aw et al.l 

l|l995l ) reveals a compact structure, while spectra and po- 
larisation measurements demonstrate that this object has 
a large UV e xcess combined with a low UV polarisation 
CPuv < 2.7 %:lTadhunter et al.ll2002l ). Moreover. [Shaw et all 
(| 19951 ') and iTadhunter et al.l (|2002l ') report the detection of 



high order Balmer lines in absorption, consistent with the 
presence of a significant YSP component. 

Spectroscopically, the classification of the AGN compo- 
nent in PKS 2135-20 is uncertain. Originally classified as an 
NLRG on the basis of its rich spect rum of narrow emissio n 
lines with moderate ionisation (e.g. ITadhunter et al.|[l993l ). 
deep EFOSC spectra reveal evidence for broad permitted 
lines (e. g. Mgll, RP) sug gesting PKS 2135-20 is instead a 
BLRG ijShaw et al.lll995h . Further, the relatively low UV 
polarisation of the source suggests that we are detecting di- 
rect rather than scattered AGN light. Similarly, the classi- 
fication of PKS 2 135-20 as a radio galaxy has been brought 
into question by ITadhunter et al.] (|l998f ) who suggest that 
its absolute magnitude of My < -23.0 is consistent with the 
classification of PKS 2135-20 as a QSO using the criteria of 
IVeron-Cettv fc VeronI (|l993l ). The presence of both quasar 
and stellar spectral features in PKS 2135-20 suggests that 
PKS 2135- 20 is a radio-loud exam ple of a post-starburst 
quasar (e.g. iBrotherton et al.lll999l ). 

Our new VLT spectrum of PKS 2135-20 is consistent 
with the detection of broad permitted lines - see Figure [8] 
for a detailed look at the features in the UV. We there- 
fore confirm the classification of this object as a BLRG or 
quasar. Given the detection of the broad lines it is likely that 
a significant power-law component is required. However, it 
is not possible to model the SED solely with an OSP and 
a power- law (x?cd ~2). Hence, we have also attempted to 
model the SED including a YSP component, consistent with 
the presence of the higher order Balmer lines in absorption 
and a significant Balmer break (see Figure |4}. In addition, 
because PKS 2135-20 has one of the highest redshifts in the 
sample, and we observe down to ~ 2300A in the rest frame, 
we have also investigated the effects of different reddening 
laws on the YSP. 

Figure [3] shows the Xrod space for three component 
(OSP, YSP and power-law) models. Consistent with the re- 
sults for the other sources studied in this paper, a three 
component model provides a large range of viable fits to the 
SED. For PKS 2135-20, a clear minimum (xLd ~ 0.45) is 
observed for YSPs with age 0.2-0.3 Gyr and reddening 
E(B-V) ^ 0.1-0.5. 

Due to the strong power-law-like continuum shape and 
broad emission lines, it is difficult to distinguish between the 
models using the stellar absorption lines - the Call K and 
the G-band are relatively weak. A few of the higher order 
Balmer lines are visible around 3800A but these appear to 
be fitted well by all of the models. A possible discrimina- 
tor could be the continuum shape between ~ 3800A and 
~ 4100A - the older age YSPs (> 0.7 Gyr) increasingly 
underpredict this region and the Call K line, whereas the 
younger age YSPs (0.04-0.6 Gyr) reproduce this region and 
Call K well. For this range of ages of YSP, the models typi- 
cally comprise: OSP (12.5 Gyr; 0-24%), YSP (0.04-0.6 Gyr, 
E(B-V) - 0.2-0.5; 45-77%) and power-law (a ~ -1.3 - -1.4; 
3-42%) with x?cd ~ 0.4-0.6. Note, the OSP and YSP con- 
tributions decrease with increasing age of the YSP, but the 
power-law component increases. An example of a model in 
the best fitting range is shown in Figure 12.5 Gyr OSP 
(22%) plus 0.2 Gyr YSP (E(B-V) = 0.2; 45%) plus power- 
law (q = -1.47; 34%) with x?ed = 0.42. 

As for some of the other sources at high redshift in our 
sample, we have also attempted to model the SED using a 
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Figure 4. SEDs and detailed fits continued 



younger OPS with age 7 Gyr — consistent with the age of 
the Universe at the redshift of the source. This gave similar 
results to the 12.5 Gyr OSP models - the level, pattern 
and location of the minimum of the x?cd space contours are 
similar, as are the detailed fits. 



4.3.12 NGC 612 

NGC 612 is the nearest radio source in our sample [z = 
0.030) and, unusually for a powerful radio source, is opti- 
cally identified with a peculiar lenticu lar/SO galaxy with a 
prom i nent dust lane an d stellar disk (jWesterlund fc SmithI 
1 19661 : lEkers et all 1 197^ . NGC 612 also contains a galaxy- 
scale disk of neutral hydrogen (H I) gas, which follows the op- 
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NGC 612: (A): 12.5 Gyr OSP + 0.1 Gyr YSP (E(B-V) = 0.9) 
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NGC 612: (B): 12.5 Gyr OSP + 0.1 Gyr YSP (E(B-V) = 0.6) 
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Figure 4. SEDs and detailed fits continued 



tical disk and dust-lane of the galaxy lEmontdliooel : Emonts 
et al. 2007 in prep). Several authors have noted similarities 
to Ce ntaurus A fe.g. lWesterlund fc Smithlll966l : lEkers et all 
Il978l ). The double radio source spans ^ 14 arcmin (~ 0.4 
Mpc) and is oriented approximately E-W, perpendicular to 
the dust lane. Due to its proximity, NGC 612 is well resolved 
spatially. Line and continuum emission is observed to large 



radii, although the emission lines (e.g. [Q II], [O III], [N II] 
and Ha) are weak (W cstcrlund fc SmithI Il966l : IGoss et"aLl 
Il980l : iTadhunter et al.l 119931 ). In addition, as well as the 
usual absorption features associated with a n OSP, the higher 
B aim er lines are detected in absorption l|Tadhunter et aU 
1 19931 . se e also Figure [H , suggesting a significant YSP com- 
ponent. iRaimann et al.l (,2005. ) have recently modelled the 
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NGC 612: (C): 12.5 Gyr OSP + 0.1 Gyr YSP (E(B-V) = 0.4) 
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NGC 612: (D): 12.5 Gyr OSP + 0.05 Gyr YSP (E(B-V) = 0.6) 
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NGC 612: (E): 12.5 Gyr OSP + 0.05 Gyr YSP (E(B-V) = 0.4) 
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Figure 4. SEDs and detailed fits continued 



stellar populations in NGC 612 using five distinct compo- 
nents, finding evidence that > 30% of the light in the nuclear 
and extended apertures is attributable to young stars with 
age 0.003-1 Gyr with the IGyr component being the domi- 
nant YSP. 

The slit position chosen for the current study was 
aligned parallel to the large-scale dust lane and stellar disk. 



Figure [T] demonstrates that the UV excess extends out to a 
radius of at least 25 arcsec from the nucleus in both di- 
rections. Interestingly, the UV excess increases with dis- 
tance from the nucleus - the nucleus itself has the largest 
1)4000 ratio. This is in contrast to the other spatially re- 
solved sources in the sample which generally have either a 
constant D4000 ratio across the galaxy or a peak in the UV 
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Table 3. Best fitting models to the galaxy continua. Columns are: (a) aperture (see Figure [T] for sizes and positions; (6) percentage 
contribution of the subtracted nebular continuum by flux in the bi n 3540- 3 640A (i.e. corresponding to the region of the maximum nebular 
continuum contribution) reddened (where applicable) using the ISeatonI l|l979l ') reddening law; (c) percentage of OSP (12.5 Gyr); (d) 
percentage of YSP; (e) age of YSP in Gyr; (/) degree of reddening required for the YSP; (g) percentage of power law component; (h) 
spectral index of power law; (i) percentage of total stellar mass attributed to the YSP component; (jr) best reduced chi squared obtained. 
* denotes galaxies best modelled with a different age of OSP: PKS 0039-44 (7 Gyr), PKS 0409-75 (7 Gyr). 

f all ages of YSP give an acceptible fit (x^cd ^'^'^ absorption lines) but there exists a clearly defined minimum in the Xrcd ^P^ce - see text. 
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excess in the nucleus. In addition to modelling the nuclear 
aperture, we have also extracted 6 extended bins in the disk 
of the galaxy. All apertures are highlighted in Figure [T] 

The nucleus. It is not possible to model the nuclear 
aperture of NGC 612 with an OSP, with or without an 
power-law component (x^d > 2) and so a YSP component is 
essential to model the SED. Figure[3]shows the Xicd space for 



two component (OSP plus YSP) models , in which the YSP 



Calzetti et al 



is reddened using both the lSeatonI l| 19791 ) and 

(200d) reddening laws. We find tha t use of 

( i200(]D rather than the ISeatonl l|l979l ) reddening law for the 



Calzetti et al 



YSP makes little difference to our results (see Section 15.2 
therefore for the analysis pres ented h e re, w e shall focus on 
the results obtained using the lSeatonI (|l979l ) law. 
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Two component fits (OSP plus reddened YSP) provide 
viable fits for ages 0.1-2.0 Gyr with two distinct minima at 
~ 0.1 Gyr (E(B-V) ~ 1.2) and ~ 0.9 Gyr (E(B-V) ~ 0.6) 
with x?ed = 0-29 and 0.31 respectively (Figure |3]). Hence, as 
with many of the apertures, it is impossible to distinguish 
between the models on the basis of the Xrod to the overall 
SED alone. 

Based on the detailed fits, models including YSPs with 
ages 0.03-0.1 Gyr (E(B-V) ~ 1.4-1.2) provide the best fit to 
the absorption features and continuum in the range 3700- 
4300A, with contributions of 80-75% and 20-28% for the 12.5 
Gyr and YSP components respectively. Models including 
YSPs with ages > 0.2 Gyr tend to either over-estimate the 
continuum in the region 4000-4100A and/or over-estimate 
the strength of the Ca K line, especially as the age of the 
YSP increases. Figure 3] shows the SED and Balmer line fits 
for the 12.5 OSP plus 0.05 Gyr YSP with E(B-V) = 1.2 and 
X?od = 0.3. 

The fits were not significantly improved when including 
a power- law component. This is not surprising given that 
NGC 612 has no strong evidence for a power law component 
- the emission lines are weak and previous work has not 
disc overed significant p olarised light (optical polarisation < 
2%: lBrindle et al.lll990^ . 

Extended apertures As discussed above, NGC 612 
is spatially well resolved and the UV excess covers the en- 
tire spatial exent of the galaxy. Further, given that an AGN 
component is not likely to be significant in the nuclear re- 
gions, it is also logical to assume the extended regions also 
do not have a power-law component. The analysis therefore 
focusses on modelling the extended aperture SEDs using a 
combination of OSP and reddened YSP components. 

As in the nuclear aperture, a pure OSP fails to provide 
an adequate fit to the SED in all of the extended apertures 
(Xred > 6) and a YSP component is required. Indeed, in- 
cluding a YSP component significantly improves the fits in 
all extended apertures - solutions are found with x?od < 1j 
see Figure [3] 

The minimum chi-squred plots for apertures A, B, D 
and E show two minima at ~0.05 - 0.1 Gyr (E(B-V) ~0.5 - 
1.0) and ~1.0 Gyr (E(B-V) ~0.0 - 0.3). However, consistent 
with fits to the spectra of other objects, the older ages (sec- 
ond minimum) can be ruled out due to the over prediction 
of the strength of CaH K, and failure to provide an adequate 
continuum fit in the region around 3700-4300A. The cases of 
apertures C and F are more clear-cut in the sense that only 
models with relatively young ages (0.04 - 0.1 Gyr) and mod- 
erate reddening (0.1 < E[B — V) < 0.2) provide adequate 
fits to the SEDs. It is notable, however, that although the 
YSP models with ages ~0.04 - 0.1 provide a good fit to the 
continuum SEDs and the higher Balmer absorption lines, 
they significantly underpredict the strength of the Call K 
absorption feature in all the extended apertures. One plau- 
sible explanation for this is that Call K absorption due to 
intervening ISM is important for these apertures. Such an 
explanation is consistent with the relatively high reddening 
deduced for the YSP in these extended apertures, as well as 
the alignment of the slit with the large-scale dust lane. 

It is also worth noting that, in apertures C and F, the 
absorption line profiles are complex and appear to have a 
split profile. At these positions, there is an enhancement in 
the continuum brightness on the 2D spectra and, in the case 



of aperture C, a 'blob' in the galaxy halo. This may be the 
signature of stellar populations which are offset in velocity 
e.g. a young star-forming cluster in the halo which is moving 
differently to the older stars in the galaxy. However, the 
data presented here have insufficient resolution and S/N to 
investigate this issue and we have assumed the absorption 
profiles can be modelled by up to two stellar populations at 
the same velocity. 

Summary for NGC 612. The high spatial resolution 
of the continuum structure in NGC 612 allows us to note 
some interesting trends in the stellar populations in this ra- 
dio galaxy: the model fits for all apertures are consistent 
with a combination of an OSP and a reddened YSP of age 
~0.04 - 0.1 Gyr, with the contribution of the YSP decreasing 
towards the nucleus — as expected for an SO galaxy. How- 
ever, although we cannot absolutely rule out a contribution 
by an older (~1 Gyr) reddened YSP in the nuclear regions, 
our spectra are not consistent with the prese nce of the dom- 
inant ^ IGyr unreddened YSP deduced bv lRaimann et ahl 
(|2005l ). 



5 DISCUSSION - DETERMINING THE 
PROPERTIES OF THE YSP IN RADIO 
GALAXIES 

The results presented in this paper provide a good illustra- 
tion of the challenges faced when determining the propo- 
erties of the stellar populations in AGN host galaxies. The 
issues that need to be considered are discussed in the fol- 
lowing sections. 

5.1 AGN-related components. 

As st ressed in previous studies (e.g. iTadhunter et al.|[l996l . 
I2002D . AGN continuum and emission line components have 
a major effect on the the continuum spectra of radio galax- 
ies at optical and UV wavelengths. In the current study, 
scattered AGN light is important for at least two objects 
(PKS 0039-44, 3C 321), a further three objects show good 
evidence for directly transmitted AGN light (PKS 1932-46, 
PKS 2135-20, 3C 381), and nebular continuum and emission 
line infilling of absorption lines is important in most objects 
(see Table |4]). To illustrate this point Figure |9] shows an 
example of the correction for both nebular continuum and 
the emission line infilling of the higher order Balmer lines in 
the case of the extended aperture of 3C 321. Not only does 
the nebular continuum make a large contribution to the UV 
continuum below the Balmer break in this case (~25%), 
but infilling of the absorption lines is also important: the 
higher order Balmer lines (e.g. H9, HIO) only become vis- 
ible in absorption following the subtraction of our nebular 
model (based on H/3). Clearly, failure to take into account 
the AGN-related components may lead to the wrong solu- 
tion in terms of modelling the YSP. 

5.2 Reddening of the YSP. 



We often find that YSP models with significant redden- 
ing are required in order to obtain good fits to the spectra of 
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Figure 9. The effect of the nebular correction in the extended 
aperture of 3C 321. The different spectra are labelled in the 
plot and the dotted lines are to aid the eye for comparisons be- 
tween the original spectrum, the nebular model and the nebular 
model subtracted spectrum. In this aperture, the nebular compo- 
nent contributes ~25% of the flux in the region 3540-3640A and 
~12% of the flux in the H9 window (3822-3858A) defined by 
iRaimann et all Hooi). The location of the H9 window and the 
WLB window (3810-3988A), both defin ed by iRaimann et al.l 
l|2005l ') and the position of HIO (used by lAretxaga et al.ll200lF 
are marked on the plot. 



many of the objects/apertures in our sample. Although the 
importance of reddening fo r the YSP h as been discounted in 
some previous studies (e.g. IRaimann et al...2005. ). the detec- 
tion of significant reddening in the current work is not sur- 
prising given the morphological evidence for dusty circum- 
nuclear environment s from HST images of several nearby 
radio galaxies (e.g . Ide Koff et all l200d : lAllen et al.1 12OO2I : 
iMadrid et al.l I2OO6I ) . It is also notable that, while the red- 
dening is often most important in the circumnuclear regions, 
it can also be significant in the extended haloes and disks of 
the galaxies (e.g. NGC 612). 

Figure 1101 hi ghlig hts the differences in shape between 
the lSeatonI \l97^ ) and lCalzetti et al.l dioool) redd ening laws. 
At optical wavelengths, the lCalzetti et aP ()200Cl ) cur ve rises 
more steeply towards longer wavelengths than the ISeatonI 
l|l979l ) law. This difference in shape causes the changes in 
the contour plots for NGC 612 (F igure[3l) - a steeper redden- 
ing law (e.g. ICalzetti et al.|[2000l ) will require l ess reddening 
in terms of E(B-V) than a shallower curve (e.g. lSeatonlllQrd l 
to model the SED with the same age YSP. The difference in 
the optical wavelength range is subtle and so does not pose 
a major problem for our results as demonstrated by Figure 
[3l However, in the UV, the c urves are mar kedly different 
with a promi nent hump in t he SeatonI (Il979l ) curve not ob- 
served in the lCalzetti et all (|2000l ) curve. Hence, the choice 
of reddening curve requires more careful consideration for 
spectroscopic studies of higher redshift objects that sample 
the rest-frame UV. 



5.3 Uniqueness of solution. 

Perhaps the greatest problem faced when modelling the 
YSP in radio galaxies is the degeneracy in the solutions 




Figure 10. Extinction law comparisons. The unreddened 0.01 
Gyr YSP from BC0 3 is plo t ted (t op) and is reddened by E(B- 
V) = 0.5 using the ISeatonI lll979h extinction law (middle) and 
the Calzct ti et al.l l|2000l ) extinction law (bottom). The extinction 
laws diffe r in two respec ts. First, the re-normalisation fact or, R {R 
= 3.1 for lSeatonI |l97i), R = 4.1 for lCalzetti et al.l |200o|)) which 
affects the overall level of the s l ope a nd second, the shape in the 
region 2OOO-25OOA fthe lSeatonl l|l979h has a prominent 'hump' at 
2175A, a feature observed in the Galaxy but rarely in extragalac- 
tic sourc es, see Section 15. 2"|l w hich affects the SED slope in the 
UV. The ICalzetti et al.l ||2OT0| ') curve also has an overall steeper 
slope towards the red. Hence, it is important to also consider the 
effect of using a particular reddening law will have on the results 
when observing sources into the UV. 



(see also iTadhunter et al.l[2005h . Despite the fact that we 
have adopted the pragmatic approach of modelling the spec- 
tra with the minimum number of components necessary 
(OSP, YSP, power-law), in some cases a wide range of YSP 
age/reddening combinations can give acceptable fits to the 
SEDs and absorption features. Unsurprisingly, the degener- 
acy problem is at its worst for objects in which the YSP com- 
ponent makes a small contribution to the optical/UV contin- 
uum relative to the OSP and/or AGN components; the most 
clear-cut solutions are obtained for the objects/apertures in 
which the YSP contributes >40% of the total flux in the 
normalising bin. While it is possible to remove some of this 
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degeneracy through detailed examination of model fits to ab- 
sorption lines that are not affected by emission line infilling 
(e.g. Call K, G-band, Mgib and, in some case, higher order 
Balmer lines), in a number of cases it has proved necessary 
to use auxilary information (e.g. UV polarisation, HST point 
source, power-law slope) to distinguish between model solu- 
tions. Such information is particularly valuable when distin- 
guishing between a very young YSP that is highly reddened, 
and a scattered or direct AGN component. 

It is also important to stress that, because it is not 
possible to completely eliminate systematic errors in the flux 
calibration, and the exact choice of bins will also affect the 
precise x?cd values of fits, the best solutions numerically (in 
terms of x?cd or other numerical measures of goodness of 
fit), are not necessarily the best solutions in terms of the 
detailed fits to continuum and absorption lines. 

Overall, it is clear that considerable care is required in 
order to avoid degeneracies in the YSP solutions; the greater 
the number of stellar components included in the fits the 
more serious the degeneracy problem is likely to become, 
especially considering that the different YSP components 
may be reddened by different amounts. This is why we pre- 
fer to model our data with the minimum number of stellar 
components, rather than assume several components at the 
outset. 



6 SUMMARY OF THE MAIN RESULTS 

Despite the challenges inherent in studies of this type, we 
have been able to extract useful information about the de- 
tailed properties of the YSP in nine out of the twelve ra- 
dio galaxies in our sample. For the remaining three galaxies 
(PKS 0039-44, PKS 1932-46 and 3C 381) we failed to find 
strong evidence for a YSP based on our spectroscopy, al- 
though this does not rule out a contribution from a YSP 
at some level. Indeed, because of the dominance of the scat- 
tered AGN component in PKS 0039-44, and the direct AGN 
component in PKS 1932-46 and 3C 381, it would be difficult 
to detect a YSP that contributes as much as ~20% of the 
total light, and a greater proportion of the stellar light, in 
such objects. We now summarise the main results for the 
YSP in the remaining ten objects for which we have been 
able to extract useful information. 



6.1 Spatial extent of the UV excess 

A striking aspect of this study is that, for all of the ob- 
jects that are well-resolved, the UV excesses are spatially 
extended, covering the full measurable extents of the galax- 
ies. Similar r esults were fo und for the four g alaxies studied 
bv [Tadhu ntc r et al.l (|2005l ) fc lEmonts et al.l ([2006). More- 
over, in all cases in which we have sufficiently high S/N to 
model the spectra of the off-nuclear emission in detail we 
find clear evidence for young stellar populations. It is also 
notable that for the two double nucleus systems in our sam- 
ple, both of the nuclei in each system show a UV excess. 
Together, these results suggest that where starbursts occur 
in radio galaxies, they tend to be galaxy-wide rather than 
confined to the immediate vicinities of the nuclei hosting the 
AGN. 



6.2 Ages of the YSP 

The ages of the YSP have the potential to provide key infor- 
mation for understanding the order of events which tr i ggere d 
the AGN in galaxy mergers (see iTadhunter et all 120051 ). 
With regard to the YSP age determination we can divide 
the current sample into four categories. 

(i) Objects with indeterminate YSP ages: 3C 236, 
PKS 2135-20, 3C 285. In these cases the contribution 
of the YSP is relatively minor over the wavelength range 
covered by the observations, and, although we believe that 
there is good evidence for YSPs in these cases, it has not 
proved possible to determine the YSP ages with any accu- 
racy: all YSP ages ranging from moderately young (~0.04 
- O.f Gyr) to intermediate (~0.5 - I.O Gyr) are possible in 
these objects. 

(ii) Objects with intermediate YSP ages (0.3 — 
3.0 Gyr): 3C 3 21. In contrast to some recent studies (e.g. 
Tadhunter et all l2005l : lEmonts etall l2006l : iRaimann et all 
2005h we find unequivocal evidence for intermediate age 
YSPs in only one object: 3C 321. In that case the YSP in 
the extended apertures (which makes a relatively large con- 
tribution to the total optical flux: 40-70%) have age ~0.7- 
f .0 Gyr. Although the model results for the other apertures 
in 3C 321 are less clear-cut, fits to the spectrum of the main 
(SE) nucleus are also consistent with an intermediate YSP 
age. 

(iii) Objects with moderately young YSP (0.05 — 
0.2 Gyr): 3C 218, NGC 612, 3C 433. Because of the 
relatively large proportional contribution of the YSP, weak 
AGN, and clearly detected Balmer absorption lines, it has 
proved possible to pin down the YSP ages accurately in 
3C 218 and NGC 612. In both of these cases the YSP are 
moderately young. 

(iv) Objects with very young YSP (<0.05 Gyr): 
PKS 0023-26, PKS 0409-75. A new finding from this 
study is the existence of a group of powerful radio galax- 
ies with very young, reddened YSP. The clearest case from 
the SED fits is PKS 0023-26, in which we detect the higher 
Balmer lines directly in absorption, but, based on consider- 
ation of auxiliary data, PKS 0409-75 is also a strong can- 
didate for an object containing a very young YSP. It is not 
surprising that the young stellar populations in these cases 
are also heavily reddened, since starbursts are expected to 
be heavily dust enshrouded in the early stages. 

6.3 Reddening of the YSP and broad-band 
colours 

A large fraction (^^70%) of the objects for which we have 
good information on the YSP show evidence for significant 
reddening in their nuclear regions (E(B-V) > 0.4), and in 
one case - NGC 612 - we also find evidence for substantial 
reddening in the off-nuclear regions. Since many studies of 
AGN host galaxies are based on the broad-band colours, it 
is interesting to consider the impact that such reddening has 
on the broad-band continuum colours of the objects in our 
sample. In Table |4] we show the (B-V) colours derived from 
our spectra for the objects in our sample along with the three 
objects in Tadhunter et al. (2005). These colours have been 
derived by measuring the ratio of the fluxes measured over 
500A bins at the central wavelengths of the B and V filters 
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Table 4. B — V colours of the nuclei in the rest frame, after 
correcting for Galactic extinction. To avoid contamination from 
the emission lines, the colours were measured using the best fitting 
models to the SED. For comparison, we include the B—V colour o f 
a 12.5 Gyr stellar population from the lBruzual fc Charloti (|2003| ) 
models. 



Object 


B-V 


Object 


B-V 


3C 218 


0.67 


PKS 0023-26 


0.93 


3C 236 


1.05 


PKS 0039-44 


0.86 


3C 285 


0.75 


PKS 0409-75 


0.84 


3C 321 


0.71 


PKS 1932-46 


0.88 


3C 381 


1.00 


PKS 2135-20 


0.52 


3C 433 


0.88 


NGC 612 


1.03 


12.5 Gyr 


0.92 







in the rest frames of the objects. For reference, quiescent 
elliptical galaxies have B — V ^ 0.9-1.0. 

From Table |4] it is clear that some galaxies have E(B- 
V) colours that are significantly bluer than quiescent ellipti- 
cal galaxies, as expected in the case of major contributions 
from unreddened YSP and/or scattered AGN components. 
However, many of the sources have broad-band colours that 
are similar to, or significant redder than, quiescent elliptical 
galaxies, despite the clear detection of UV-excesses based 
on 4000A, break measurements and, in some cases, strong 
evidence for YSPs from spectral fitting and the detection of 
Balmer absorption lines. These results serve to demonstrate 
that considerable caution is required when making conclu- 
sions about the stellar populations in AGN host galaxies on 
the basis of broad-band colours alone. 



6.4 The mass contribution of the YSP 

Another significant feature of our results is that the YSP 
typically make up a significant fraction of the total stellar 
mass in the regions sampled by the spectroscopic slits (~1 
- 35%). Note that this is true even for some of the objects 
in which the measured age of the YSP is relatively young. 
In such cases the YSP are often heavily reddened, therefore, 
although the YSP contribute a relatively small fraction of 
the total hght, they can represent a significant fraction of 
the total stellar mass. 



7 CONCLUSIONS 

We have modelled deep wide spectral coverage data for a 
sample of 12 powerful radio galaxies with previous evidence 
for young stellar populations. The results demonstrate the 
importance of accounting for AGN contamination and red- 
dening when determining the detailed properties of the stel- 
lar populations. Moreover, in cases where the YSP makes a 
relatively small contribution to the optical/UV continuum, 
it has proved difficult to determine the precise ages of the 
YSP. We find that the YSP with well-determined properties 
cover a wide range of post-starburst ages (0.03 - 1.5 Gyr) 
and reddening (0.0 < E{B — V) < 1.5), and account for 
a significant fraction of the total stellar mass in the regions 



covered by our spectra (1 — 50%). The implications of these 
results for our understanding of the evolutionary status of 
radio source host galaxies, and the heating mechanism of the 
warm/cool dust in AGN, will be discussed in forthcoming 
papers (Tadhunter et al. 2007a,b). 
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